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Abstract. In this paper we present a general procedure that allows for the 
reduction or expansion of any network (considered as a weighted graph) . This 
procedure maintains the spectrum of the network's adjacency matrix up to a 
set of eigenvalues known beforehand from its graph structure. This procedure 
can be used to establish new equivalence relations on the class of all weighted 
graphs (networks) where two graphs are equivalent if they can be reduced to 
the same graph. Additionally, dynamical networks (or any finite dimensional, 
discrete time dynamical system) can be analyzed using isospectral transfor- 
mations. By so doing we obtain stronger results regarding the global stability 
(strong synchronization) of dynamical networks when compared to other stan- 
dard methods. 

1. Introduction 

Real world networks, i.e. those found in nature and technology, typically have 
a complicated irregular structure and consist of a large number of interconnected 
dynamical units. Coupled biological and chemical systems, neural networks, social 
interacting species, and the Internet are only a few such examples [51 UM [T51 [Ml 121)1 
[28] . Because of this complexity, the first approach to capturing the global properties 
of such systems has been to model them as graphs whose nodes represent elements 
of the network and the edges define a topology (graph of interactions) between these 
elements. This is principally a static approach to modeling networks in which only 
the structure of the network's interactions (and perhaps the interaction's strength) 
are analyzed. 

A more general line of research is concerned with the dynamical properties of 
such networks. Mathematically, this has been done by modeling networks as inter- 
acting dynamical systems [TJ HI [T2] . Important processes that are studied within 
this framework include synchronization or contact processes, such as opinion 
formation and epidemic spreading. These studies give strong evidence that the 
structure of a network can have a substantial impact on the network dynamics [5] . 

Regarding this connection, the impact of a dynamical network's underlying struc- 
ture on its dynamics is of major interest. In this the spectrum of the network's 
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adjacency matrix has naturally emerged as a key quantity in the study of a variety 
of dynamical networks. For example, systems of interacting dynamical units are 
known to synchronize depending only on the dynamics of the uncoupled dynamical 
systems and the spectral radius of the network's adjacency matrix [20j [2T] . More- 
over, the eigenvalues of a network are important for determining if the dynamics of 
a network are regular or chaotic. [U H] . 

To aid in understanding this interplay between the structure (graph of inter- 
actions) and dynamics of a network in this paper we introduce the concept of an 
isospectral graph transformation. Such transformations allow one to modify the 
structure of a network, at the level of a graph, while maintaining properties related 
to the network's dynamics. For the sake of the reader, we will present several of 
our results regarding isospectral graph transformations here in the introduction. 
However, we will do so in a very informal manner as exact formulations require the 
introduction of many formal definitions. 

Mathematically, an isospectral graph transformation is a graph operation that 
modifies the structure of a graph while preserving the eigenvalues of the graph's 
(weighted) adjacency matrix. Besides modifying interactions, such transforms can 
also reduce (or increase) the number of nodes in a network. As not to violate the 
fundamental theorem of algebra, isospectral graph transformations preserve the 
spectrum of the graph (in particular the number of eigenvalues) by allowing edges 
to be weighted by rational functions. In other words, a matrix of smaller size can 
have the same spectrum as a larger matrix if the entries of the smaller are allowed 
to be functions of a spectral parameter. This fundamental result is contained in 
theorem 13.51 

For a graph G with vertex set V this is done by reducing (or expanding) G with 
respect to specific subsets of V known as structural sets (see section 3.1 for exact 
definitions). As a typical graph has many different structural sets it is possible 
to consider different isospectral transformations of the same graph G as well as 
sequences of such transformations. However, given the large number of possible 
reductions of a typical graph the results of theorems 13.81 and 13.91 are quite useful 
and central to the theory of isospectral graph reductions. 

Theorem 3.8 and 3.9 (Existence and Uniqueness). If G is a graph with vertex 
set V then it is possible to sequentially reduce G to a graph on any nonempty subset 
V C V of its original vertex set. Moreover, the resulting graph does not depend on 
the number or sequence of intermediate reductions but only on the choice of the 
set V. 

Because it is possible to uniquely reduce a graph over a subset of its vertices, 
isospectral graph reductions can be used to induce new equivalence relations on the 
set of all graphs. 

Theorem 3.11 (Spectral Equivalence). Suppose we have a rule r that selects 
a unique subset of vertices t(G) C V for any graph G. Then the rule r induces an 
equivalence relation (spectral equivalence) on the set of all graphs (or all networks 
considered as graphs). Namely, two graphs (networks) are spectrally equivalent if 
they can be isospectrally reduced to the same graph (network). 

Observe that in order to induce an equivalence relation we need a rule that selects 
a unique subset of the vertices of a network. For instance, the rule "choose any 
vertex" does not work. Given this restriction, an expert in biology, sociology, etc. 
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must determine appropriate vertex sets or rules over which to reduce e.g. minimal 
(maximal) degree, minimal (maximal) clustering coefficient, etc. Once this is done 
our procedure can then be used to reduce the network(s). However, the rule must 
be chosen by an expert in a corresponding applied field. 

Importantly, if using a particular rule does or does not yield insight into a par- 
ticular network (or set of networks) it is possible to try another rule. In this regard, 
our procedure provides a working tool for scientists, engineers, and others for the 
analysis of their respective networks. However, determining the correct rule to use 
requires an experts knowledge of his/her own field. 

Network reductions themselves are motivated by the need to find ways of re- 
ducing a network's complexity while simultaneously maintaining some important 
network characteristic(s). This amounts to coarse graining or finding the right scale 
at which to view the network. Isospectral network transformations, specifically re- 
ductions, serve this need by allowing one to view a network as a smaller network 
with essentially the same spectrum. Finding the right scale here then amounts to 
finding the right structural set or rule over which to reduce the network, which is 
the job of the expert. 

We note that by allowing the edge weights of the graphs to be functions it may 
appear that our procedure is trading the complexity of the graph's structure for 
complex edge weights. In fact, it is often possible to isospectrally transform a graph 
in a way that does not change edge weights of the graph or preserves the type of 
edge weights a graph has e.g. positive integers, real numbers, etc. (See theorems 
Oandlil) 

Moreover, it is worth noting that one can reduce the line graph (often called an 
edge graph) L(G) of a graph G. As the line graph of L{G) is formed by associating 
a vertex with each edge of G then reducing over a structural set of L{G) amounts 
to reducing G over a particular edge set. From this point of view, isospectral graph 
reductions can be done either over vertex sets or edge sets. 

Of primary interest is the fact that isospectral graph transformations suggest 
other useful transformations on interacting dynamical systems i.e. dynamical net- 
works. An important example in this paper is a dynamical network expansion in 
which a dynamical network is modified in a way that preserves its dynamics but 
alters its associated graph structure. 

Such transforms provide a new tool for the study of the interplay between the 
structure (topology) and dynamics of dynamical networks. Much as general dy- 
namical systems are investigated via change of coordinates, isospectral network 
transforms introduce a mechanism for rearranging the specific (graph) structure of 
a system while preserving the networks dynamics in an essential way. By so doing 
the original network's dynamics can be investigated by studying the transformed 
dynamical network. The main result of this technique is the following. 

Theorem 6.6 and 6.8 (Dynamical Network Expansions). If an expansion of 
a dynamical network has a globally attracting fixed point then the network itself 
has a globally attracting fixed point. Moreover, determining whether a dynamical 
network has a global attractor is more easily done by analyzing an expansion rather 
than the original network. 

That is, the method of dynamical network expansion provides a new efficient tool 
for establishing the global stability of dynamical networks. Moroever, the result(s) 
theorem [6. 5 1 and theorem [6.71 generalize the results given in [HE]. Moreover, as any 
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finite dimension dynamical system with discrete time can be considered a dynamical 
network this technique can be applied to a very large class of dynamical systems. 

The structure of the paper is as follows. In section 2 we introduce dynamical 
networks along with their graph of interactions. We then give sufficient conditions 
under which a dynamical network will have a globally attracting fixed point. Sec- 
tion 3 introduces the procedure of an isospectral graph reduction and results on 
sequences of such transformations. Graph transformations that preserve weight 
sets and spectral equivalence are treated in section 4. The proof of the theorems 
contained in sections 3 and 4 are then given in section 5. In section 6 the behavior 
of dynamical networks is investigated via dynamical network expansions. Section 
7 contains some concluding remarks. 



2. Preliminaries 

Dynamical networks or networks of interacting dynamical systems are composed 
of (i) local dynamical systems which have their own (local intrinsic) dynamics, (ii) 
interactions between these (elements of the network) local systems, and (iii) the 
graph of interactions (topology of the network) . 

2.1. Dynamical Networks and their Graph Structure. Our first task is to 
establish a mathematical framework for the investigation of dynamical networks. 
This is done following the approach in [1] . 

Let i G I = {1, . . . , n} and Tj : Xi — >• Xi be maps on the complete metric space 
(Xi, d) where 

/ -i \ T d(T l (x l ),T l (y l )) 

(1) Li = SUp j- r < 00. 

xijtyiGXi a(x l ,y i ) 

Let (T, X) denote the direct product of the local systems (T,-, JQ) over L on the 
complete metric space (X, d max ) where for xjgX 

d m ax(x,y) = max{d(a; i ,?/ i )}. 

Definition 2.1. A map F : X — > X is called an interaction if for every j € X there 
exists a nonempty collection of indices Ij C I and a continuous function 

F . :0.V, 

ieij 

that satisfies the following Lipschitz condition for constants Ay > : 

(2) d(F, J ({x t }),F J ({y l })) < Aydfo.i/i) 

for all {xi}, {yi} G Q) ieI ^ where {xi} is the restriction of x G X to © ieI . Xi. 
Then the (interaction) map F is denned as follows: 

=F j ({x i }), j el, ieij. 

The constants Ay in definition 12.11 form the Lipschitz matrix A G M" x " where 
the entry Ay = if i £ Ij . 



Definition 2.2. The superposition T = F o T generates the dynamical system 
(T, X) which is a dynamical network. 
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Figure 1. The graph Tp corresponding to the interaction F in 
example [T] 



To each dynamical network (F, A) there is a corresponding graph of interactions. 
A graph of interactions is an unweighted directed graph representing the structure 
of the interaction F between the network elements. 

An unweighted directed graph G is an ordered pair G = (V, E). The sets V and E 
are the vertex set and edge set of G respectively. If the vertex set V = {v\ 1 . . . , v n } 
then we denote the edge from v% to Vj by . 

Definition 2.3. Let F : X —> A be an interaction. The graph Tp = (V,E) with 
vertex set V = {vi, . . . ,v n } and edge set E = {eij : i G Ij, j el} is called the 
graph of interactions of F. 

We note that each vertex Vi G V of the graph of interactions Tp = (V,E) 
corresponds to the ith. element of the dynamical network (_F, A). Moreover, there 
is an edge G E if and only if the jth coordinate of the interaction -F(x) depends 
on the ith coordinate of x. 

Remark 1. Suppose the interaction F : X — > X is continuously differentiable and 
each Xj C M. If DF is the matrix of first partial derivatives of F then the constants 

Aij = max \(DF)u(x)\ 
satisfy condition © for the interaction F. 

Example 1. Let F : [0, l] 3 — > [0, l] 3 be the interaction given by 

Xi(l - ax 2 ) 
F(x) = fixixz 

_ x 3 {l - -yx 2 ) 

Using the Lipschitz constants A^ = max l£ x \ (DF)ji(x)\ the interaction F has 
Lipschitz matrix 

"1/30" 
A = a 7 
0/3 1 

The graph of interactions T p of F is the graph shown in figure [7J 



for <a,£,7 < 1. 



2.2. Stability of Dynamical Networks. One of the goals of this paper is to find 
sufficient conditions under which a dynamical network (J 7 , A) has simple dynamics. 
By simple we mean that the system (J 7 , A) has a globally attracting fixed point. 
That is, that all trajectories of the dynamical network converge to a single point as 
time tends to infinity. 

Definition 2.4. The dynamical network (F, A) has a globally attracting fixed point 
x G A if for any x G A, 

lim d ma J.F fc (x),x) = 0. 

k— ►oo 
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If (J 7 , X) has a globally attracting fixed point we say it is globally stable. 

Recall that the constants Li and Ay- come from (fTJ) and @ for the local systems 
(T, X) and interaction F respectively. For the dynamical network (J 7 , X) we define 
the matrix Mjr = A T ■ diag[Li, . . . , L n ]. That is, 



/ AuLi . . . K nl L n 



\ Ai„Li . . . A nn L r , 
Wc let p(Mjr) denote the spectral radius of the matrix A, i.e. if a(Mjr) are the 
eigenvalues of A then 

p(Af» =max{|A| : A G tr(M^)}. 

The following theorem gives a sufficient condition under which a dynamical network 
has a globally attracting fixed point. 

Theorem 2.5. J/ p(Mjr) < 1 t/ien £/ie dynamical network (J 7 , X) has a globally 
attracting fixed point. 

Before proving theorem 12.51 we observe that if v G C" xl then its £°° norm 
|v||oo = maxi \vi\. The £°° norm of a matrix A e C nx ™ is 



3=1 

or the maximum absolute row sum of A. Moreover, the £°° matrix norm is sub- 
multiplicative, i.e. IHABIIloo < 1 1 |A| | |oo 1 1 |S| | |oo for any A, B 6 C nxn . A proof of 
theorem 12.51 is the following. 

Proof. For xjel and 1 < j < n 

d^x^y),) =d(i^({T(x) i }),i ? i({T(y)i})) 

^A^T^),^*)) 



<y^ KijLid{xi,yi). 



Therefore, each entry of the column vector 

df^x)!,^)!) 



d(j 7 (x)„, J 7 (y) r 



< M 



.f 



d(xi,yi) 









' d^x^JXy)!) " 




d(x ll y 1 ) 


As 






_ d(j-(x) n ,^(y) n ) _ 




_ d(x n ,y n ) 



by the same reasoning then inductively 
(3) d max {T k ^),T k {y))< 



M k 



d(xx,yi) 
d(x ni y n ) 
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for all k > 0. 

By the Jordan canonical form theorem there is a non-singular matrix A G C" xr 
and a block-diagonal matrix J G C™ x " such that Mjr = AJA^ 1 . In particular, 



J = 



J mi (Ai) 
J m2 (A 2 ) 













J rot (A t ) 



where 



-^mi(Ai) — 



Ai 1 

A, 1 



... 



... 



A, 1 
Xi 



G C 



mi xrrii 



and Xi £ a(Mj^) for 1 < i < t. Moreover, as J is block diagonal then for k > 0, 

J* t (Ai) ... 

o j*,(a 2 ) ; 



j fe = 







... J*, (A*) 

As the norm of a matrix is its maximum absolute row sum then 
|||J fe |||oo = max |||J^(A0|||oo. 

A well-known result states that the kih power of any Jordan block J mi (Af) is 



( k \ \ fc-mi+1 
I k \ it-mi+2 







( fe )A* 

for any k > mi — 1. Hence, for k > mi — 1 the matrix norm 

mi — 1 



l^m ( (A» 



= E 





A 



fe-j 



Ay 



m.i-2 I 



mi - 1 



ifc-mj + 1 1 



fc-mi + 1 



where C m< =ESo |(7) A 



m(-l I fk\ ,mj-j-l I 



Suppose p(Aljr) < 1. Letting m — max mi and C = max C mi then 

0<i<i Ki<t 



< Cp(Mjr) k - m+1 for all fc > m - 1. 



s 
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Via equation 

oo oo 

^d mM (^ fe (x),^ fe (y)) < \\M. 

k—rn k—m 
oo 

< 2 IIUj^'A- 1 



k—m 



< 



< 



A 



A- 1 



A- 1 



d(xi,yi) 

d(x ni y n ) 
d(xi,yi) 

d{x ni y n ) 
d(xi,yi) 

d{x n ,y n ) 
d(xi,yi) 

d(x n ,y n ) 



T,\\\ jk 

k—rn 



< oo. 



Letting y = F(x) it then follows that Y^kLm dmax (F k (x), F k+1 (x)) < oo. Hence, 
for e > there exists an N such that Y1^=n dmax (.F fc (x), F k+1 (x)) < e. Therefore, 

d max (F k ( X ),F e (x)) <e 

for any k,£ > N or the sequence {F k (x)}k>i is Cauchy. Since X is complete then 
this sequence converges. 

Suppose T (x) —5- x. As J 7 is continuous it follows that .F(x) = x. More- 
over, given that YX=\ d ma x (F k (x) , T k (y)) < oo for any y E X then F k (y) -> x 
completing the proof. □ 

Example 2. Let Tj : [0,1] — > [0,1] be the logistic map Ti(xi) — 4x^(1 — xi) for 
1 < i < 3. Le£ F : [0, l] 3 — » [0, l] 3 fee the interaction given by 

1 - x 1 x 2 /9 
F(x) = 1 - 0:1X3/9 
1 - x 2 x 3 /9 

One can check that the constants Li = 4 and 



Ay = 



2/9 i/i+jV4 
otherwise 



satisfy 0) and (0j respectively. Ln particular, Ay = max xg jc \(DF)ji(x)\. For these 
constants the matrix 

' 4/9 4/9 
Mjr = 4/9 4/9 
4/9 4/9 

As p(Mjr) = 8/9 < 1 then the dynamical network (J-,X) has a globally attracting 
fixed point. 

The local systems (T, X) are said to be stable if max je i{Li} < 1 and are said 
to be unstable otherwise. We say the interaction F : X — > X stabilizes the local 
systems (T, X) if the local systems are unstable but the dynamical network (J 7 , X) 
has a globally attracting fixed point. If the local systems (T, X) are stable and 
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(T,X) has a globally attracting fixed point we say the interaction F : X — > X 
maintains the stability of (T, X). The following is a corollary to theorem 12.51 

Corollary 1. For the dynamical network (F,X) let inaxj^z {Li} = L. If Lp(A) < 1 
then the interaction F stabilizes (or maintains the stability of) the local systems 
(T,X). 

Before proving this we require the following. For A,B€ W xn we write 

A <B if Aij < Bij for 1 < i, j < n. 

If < A < B then a well known result in matrix theory states that p(A) < p(B) 
(see chapter 8, [H]). We now give a proof of corollary [TJ 

Proof. As < A T • diag[L u . . .,£„] < LA T then p(Mjr) < p(LA T ). Given that 
LA T = {LX : A £ <r(A T )} then this implies p(LA T ) = Lp(A) which in turn implies 
that p(Mjr) < Lp{A). Hence, if Lp(A) < 1 then by theorem 12.51 the interaction F 
stabilizes or maintains the stability of the local systems (T, X). □ 

For the dynamical network T = F o T in example [2] note that L — 4 and 
p(A) — 2/9. As Lp(A) < 8/9 then corollary [T] implies that the interaction F 
stabilizes the local systems (T, X). 

2.3. Dynamical Networks Without Local Dynamics. As defined in section 
2.1 a dynamical network J- = F o T is the composition of the network's local 
dynamics T and the interaction F. However, if the system has no local dynamics, 
i.e. T = id is the identity map, then the dynamical network J- is simply the 
interaction F. 

Remark 2. Note that any composition T — F o T can be considered to be an 
interaction. Writing J- — {F o T) o id the dynamical network (J 7 , X) is simply the 
interaction F o T with no local dynamics. 

The following is another corollary of theorem 12.51 

Corollary 2. Suppose the constants Li and Ajj satisfy (QJ) and 0) for the map 

J- = F oT . Let J- = (F oT) oid be the map inducing the dynamical network (F,X) 
with interaction F o T and no local dynamics. Then there exist constants Li and 
Ay satisfying (Qp and (0) for id and F o T respectively such that p(M-p) < p(Mjr). 

Proof. Let Li and Aij be as in the statement of corollary [5] Suppose x, y G X. 
Then 

d(id(x)i,id(y)i) = d(xi,yi). 
Hence, the constants Li = 1 satisfy ([TJ for the local systems (id, X). Moreover, 

d{(F o T)(x),, (F o T)(y),) =d(F j ({T i (x i )}), ^({T^)})) 

A lj L. l d(x l ,y i ). 

iex 

Therefore, the constants A^ = AijLi satisfy equation ([2j for the interaction 
F o T. For this choice of constants note that the matrix M.p is equal to the matrix 
M.jr. This implies the result of the corollary. □ 
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Example 3. Let Ti : [0,1] — > [0,1] be the map Ti(x) — sin(7TXj) for 1 < i < 2. 
Suppose F : [0, l] 2 — > [0, l] 2 is the interaction given by 

1 - .x 2 /4 
1 - .x 2 /4 



F(x) = 



Ao£e i/iai i/ie constants Li — n and Ay- = max Ig x K-D^Oji^)! satisfy (QP ano 
(0) /or T and i 7, respectively. If the dynamical network F = F o T then 



vr/2 
tt/2 

^4s p(Mjr) = 7r/2 > 1 i/ien theorem \2.5\ does not directly apply to the dynamical 
network F. 

However, suppose J- = (F o T) o id is the dynamical network with interaction 
FoT and no local dynamics. Let Li = 1 and A.y = max ie j( | (D(F oT)) (x)\ which 
satisfy ([7]) and /or T — id and FoT respectively. With this choice of constants 

tt/4 
vr/4 



^4s p(Mjr) = 7r/4 < 1 t/ien the dynamical network (J-,X) has a globally attracting 
fixed point. Importantly, given that T = T then the dynamical (T , X) must also 
have a globally attracting fixed point. 



Note that if Li = 1 for all i then A • diag[Lx, . . . , L r 
to the following remark. 



A. This observation leads 



Remark 3. If (F, X) is considered as a dynamical network with no local dynamics 
then yVfjr = A for any constants Ajj satisfying ([2]). Hence, (F, X) has a globally 
attracting fixed point if p(A) < 1. 

By considering a dynamical network in a modified form, i.e. as an interaction, it 
is therefore possible to obtain improved estimates of the network's global stability. 
The result of corollary[2]is then the first step toward answering the following general 
question: How can a dynamical network (equivalently interaction) be modified to 
get improved estimates of a network's global stability? 

Additionally, in this paper we consider whether it is possible to transform (sim- 
plify) a dynamical network at the level of its graph of interactions while maintaining 
the dynamic properties of the network. Both questions serve to motivate the graph 
transformations introduced in the following sections. 



3. ISOSPECTRAL GRAPH REDUCTIONS 

In this section we formally describe the isospectral reduction process of a graph. 
We then give specific examples of this method and present some results regarding 
sequences of isospectral reductions and spectral equivalence of graphs. 

3.1. The class of graphs G. We consider the following class of graphs, namely 
those graphs that are weighted, directed, with edge weights in the set W[A] (defined 
below). Such graphs form the class G. A graph G G G is then an ordered triple 
G = (V, E, lo) where V and E are again the vertex set and edge set of G respectively. 
The function u : E — >• W[A] gives the edge weights of G. We adopt the standard 
convention that each edge of the weighted graph G has a nonzero weight. 
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Remark 4. The class of graphs G is very general as it contains, for instance, the 
class of undirected graphs with numerical weighs. 

For notational convenience we assume that each vertex set V is given the labelling 
V = {vi, . . . ,v n }. Also, any graph G written as an ordered triple G — (V,E,ui) 
will be implicitly assumed to be a graph in G. 

Let C[A] be the set of polynomials in the complex variable A with complex 
coefficients i.e. the polynomial ring over C with indeterminate A. We denote by 
W[A] the field of fractions of C[A] or the set of rational functions of the form p/q 
where p, q G C[A], and q ^ 0. Our first task is to extend the definition of the 
eigenvalues of a matrix A E C nxn to the eigenvalues of a matrix A(X) G W[A]" X ™. 

The element a of the set A that includes multiplicities has multiplicity m if 
there are m elements of A equal to a. Suppose A and B are sets that include 
multiplicities. If a € A with multiplicity m and a G B with multiplicity n then 

(i) the union A U B is the set in which a has multiplicity m + n; and 

(ii) the difference A — B is the set in which a has multiplicity m — nifm — n > 
and multiplicity otherwise. 

If a square matrix A G £7 IX ™ then its characteristic polynomial det (A — A/) G 
C[A]. However, if A(X) G W[A] nXTl then 

(4) det (A(X) — Xl) S W[A]. 

i.e. this determinant may be a rational function of A. Despite this we will adopt 
the convention of referring to (j4j as the characteristic polynomial of A(X). 

Definition 3.1. For a matrix A(X) G W[A] nx ™ suppose the characteristic polyno- 
mial det(A(A) — XI) = p/q where p,q G C[A]. Define the sets 

P = {X G C : p = 0} and Q = {A G C : q = 0} 

where these sets includes multiplicities. We call the sets 

<r(A(X)) = P - Q and a^ 1 (A(X)) = Q-P 

the spectrum (or set of eigenvalues) of A(X) and the inverse spectrum of A(X) 
respectively. 

Of primary importance is the fact that the representation of p/q G W[A] is not 
unique. That is, p/q is equivalent to r/s for p,q,r,s G C[A] if ps = qr. It is 
therefore necessary to show that the spectrum and inverse spectrum of a matrix 
A(X) G W[A]" xn are well defined. 

That is, p/q is equivalent to r/s in W" x ™ for p,q,r,s G C[A]. Let P and Q be 
defined as in definition 13. II Similarly define 

R = {A G C : r = 0} and S = {A € C : s = 0}. 

As ps = qr it follows that PUS = QUR. Moreover, &sPUS-Q = Ql)R-Q = R 
then (P U Q) - S = R - S. 

From the fact that (PUS-Q)-S = (PUS-S)-Q it follows that P-Q = R-S. 
Similarly, one can show that Q — P = S — R. That is, the spectrum and inverse 
spectrum given by definition 13.11 are well defined. 

Before proceeding we note the following. According to definition 13.11 any matrix 
A G C" x ™ has a well defined spectrum a (A) and inverse spectrum a^ 1 (A). In 
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particular, if A € C nxn the characteristic polynomial of A is the ratio p/q G W[A] 
where p — det(A — XI) and q = 1. Hence, 

a(A) = {X G C : det(A - A/) = 0}, ct~ V) = 0. 

This can be summarized as follows. 

Remark 5. The eigenvalues of a square matrix with entries in W[A] is a gen- 
eralization of the standard notion of the eigenvalues of a complex valued square 
matrix. 

Suppose G — (V, E, lu) with vertex set V = {vi, . . . , v n }. The matrix M(G) G 
W[A]" X " defined entry wise by 

M{G) lJ =u{e l] ) 

is called the weighted adjacency matrix of G. For the graph G we denote by er(G) 
and cr _1 (G) the spectrum of and inverse spectrum of M{G) respectively. 

3.2. Isospectral Graph Reductions. Our next task is to develop the mathe- 
matical framework, i.e. definitions and notation, needed to describe an isospectral 
reduction of a graph. This requires the following standard terminology. 

A path P in the graph G — (V,E,lu) is an ordered sequence of distinct vertices 
vi,...,v m G V such that ej,j + i G E for 1 < i < m — 1. We call the vertices 
«2, ■ ■ ■ , Vm~i of P the interior vertices of P. If the vetices v\ and v m are the same 
then P is a cycle. A cycle v± . . . , v m is called a loop if m = 1. Note that as Vi, Vi is 
a loop of G if and only if en G E we may refer to the edge en as the loop. If S C V 
where V is the vertex set of a graph we will write S = V — S. 

The main idea behind an isospectral reduction of a graph G — (V, E, w) is that 
we reduce G to a smaller graph on some subset S C V. The sets S for which this 
is possible are defined as follows. 

Definition 3.2. Let G = (V,E,ui). A nonempty vertex set S C V is a structural 
set of G if 

(i) each cycle of G, that is not a loop, contains a vertex in S; and 

(ii) ui(eu) ^ A for each Vi G S. 

Observe that assumption (i) of definition 13.21 implies that a structural set S of 
G depends intrinsically on the structure of G. On the other hand, part (ii) of 
definition 13.21 is the formal assumption that the loops of the vertices in S, i.e. the 
complement of S, do not have weight equal to A G W[A]. For G G G let st(G) 
denote the set of all structural sets of the graph G. 

Definition 3.3. Suppose G = (V, E, lu) with structural set S = {vi, . . . , v m }. Let 
Bij(G; S) be the set of paths or cycles from Vi to Vj with no interior vertices in S. 
We call a path or cycle (3 G Bij(G; S) a branch of G with respect to S. We let 

B S {G)= |J B«(G;S) 

l<i,j<m 

denote the set of all branches of G with respect to S. 

If /3 = V\, . . . , v m is a branch of G with respect to S, i.e. /3 G Bs(G), and m > 2 
then let 
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g n s (G) 

Figure 2. Reduction of G over S = {vx,V3} where each edge in 
G has unit weight. 



For m = 1, 2 let V^iP) = uj(e\ m ). We call Vu(£>) the branch product of /3. Note that 
assumption (ii) in definition 13.21 implies that the branch product of any j3 G Bs(G) 
is always defined. 

In a procedure we term an isospectral graph reduction we replace the branches 
of a graph with edges. The following definition specifies the weights of these edges. 

Definition 3.4. Let G = (V,E,w) with structural set S = {v\ ...,v m }. Define 
the edge weights 

(6) M e u) = S 0eB i} (G;S) for 1 < i,j < m. 

! otherwise 

The graph lZs(G) = (S, £, /i) where G £ if ^(e^) ^ is the isospectral reduction 
of G over 5. 

Observe that ^(e^) in definition 13.41 is the weight of the edge e,j in lZs{G). 
Moreover, as W[A] is closed under both addition and multiplication then the edge 
weights fJ,(eij) of lZs(G) are also in the set W[X\. Hence, the isospectral reduction 
IZs(G) is again a graph in G. 

Example 4. Consider the graph G = (V,E, uS) given in figure^ (left) in which 
each edge of G is assumed to have unit weight. The vertex set S = {^1,^3} C V is 
a structural set of of G as 

(i) the three nonloop cycles of G, namely v%, V2, V3, U4, v%; Wi,«5,Vi/ and W3,wg,U3 
each contain a vertex in S; and 

(ii) the loop weights vertices in S — {V2, V4, V5, vq} are ui{e2i) = 1, ^(644) = 0, 
w(e55) = 1, and uj(eee) = 0. That is, io(ea) ^ X for each Vi G S. 

Conversely, the vertex set T = {vi,vs} is not a structural set of G. This follows 
from the fact that the (nonloop) cycle vs,vq,V3 does not contain a vertex of T . 

The branches of G with respect to S can be seen to be Bi±(G;S) — {vi, V5, Vi}, 
Bis(G; S) = {vi,V2,V3}, B3\(G;S) = {f3,W4,wi}, and the two branches in the set 
B 33 (G;S) = {v 3 ,v 3 ;v3,v 6 ,v 3 }. 

Using equation $5$ the branch product of each of these branches are given by 

Vu(vi,VS,Vl) = V U (V1,V2,V3) = j— j, 

Vu,(v3,V4,Vi) =V u (v 3 ,Ve,V3) = -, and V cu {v 3 ,v 3 ) = 1. 

The weight of each edge oflZs{G) = (S,£,fJ>) given by (OJ) can then be computed 
to be 

Men) = jti(ei 3 ) = -, 
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Figure 3. Restriction of the graph G to S = {t>2, «4, «5, v§}. 



v 6 , 



M e 3i) = j, and n{e 33 ) = 1 + - = —j— . 

Note that as each edge weight is nonzero then the edge set ofTZs(G) is given by 
£ = {en, ei3, e^i, 633}. The graph lZs{G) is shown in figure^ (right). 

Recall that if S is a structural set of the graph G G G then the isospectral 
reduction lZs(G) is also a graph in G. Hence, both G and lZs{G) have well-defined 
spectra. The relation between the spectrum a(G) and a(lZs{G)) is given in the 
following theorem. 

Theorem 3.5. Let S be a structural set of the graph G G G. Then 

a(K s (G)) = (a(G) U a~ l (G\ B )) - a(G\ s ). 

The proof of theorem 13.51 will be postponed until section 5. However, we note 
here that the characteristic polynomial of the matrix M(G\§) has a particularly 
simple form. 

Since the only cycles the graph G\§ contains are loops then the following holds. 
The vertices of G\§ can be ordered such that the matrix M(G\g) is triangular i.e. 
either all the entries below or all the entries above the main diagonal of M(G\ S ) 
are zero (see Frobenius normal form [S]). Hence, the matrix M(G\ S ) — XI is also 
triangular. 

As the determinant of a triangular matrix is the product of its diagonal entries 
then 

(7) det (M(G\s) - M) = J] H e «) - A )- 

Because of the simplicity of equation both the spectrum <j(G|g) and inverse 
spectrum cr^^GIg) can be calculated with relative ease. Moreover, in light of 
equation ([7]) theorem 13.51 has the following corollary. 

Corollary 3. Let S be a structural set of the graph G G G. If M{G) G C nxn then 

(i) cr(G\ s ) = {uj(eu) : Vi G 5}; 

(it) o- x (G\s) = 9; and 

(in) <r{K s (G)) = tr(G) - a(G\§). 

In applications the graphs (matrices) that are used often have real or positive 
weights (entries). If G = (V, E, u>) has complex valued weights and S G st(G) then 
corollary [3] states the following. The spectrum of TZs(G) and G differ at most by 
the spectrum of G\§. Moreover, the spectrum a(G\ s ) can be described as follows. 

For Vi G S the weight uj(eu) G C is the weight of the loop en if en is an edge of 
the graph G. In the case that en £ E then ui(eu) = 0. That is, a(G\§) consists of 
all the weights of the loops off the structural set S. This includes the weight(s) 
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for all Vi € S such that en ^ E. Because of this, the set <j(G\§) can be easily read 
off the graph G. As an illustration we give the following example. 

Example 5. Let G be the graph considered in example^ As previously shown the 
vertex set S — {vi,vs} is a structural set of G. Moreover, M(G) £ C 6x6 . Hence, 
corollary^ allows us to quickly compute the eigenvalues of the reduced graph lZg{G) 
once the eigenvalues of G are known. 

As one can calculate, the eigenvalues of the graph G are o~{G) — {2, — f , f , f , 0, 0}. 
The restriced graph G\§ shown in figure\3\has loop weights w(e22) = 1, ^(e^) = 0, 
w(e55) = 1, anduj(eee) = 0. Then corollary\3\implies that a{G\g) — {1,0, f,0} (and 
a~ 1 (G\s) — %). Additionally, as a(lZs(G)) = cr(G) — a(G\§) then the spectrum of 
the reduced graph is a(lZs{Gy\ = {2, —1}. 

That is, in the reduction process of G over S we lose the eigenvalues {0,0, 1, 1}. 
However, even if o~(G) is unknown we still know the following. The set of eigenval- 
ues o~(G\§) — {0, 0, 1, 1} is the most by which a(lZs(G)) and o~(G) can differ. 

Remark 6. We note that it is possible when reducing a graph to maintain the 
graph's entire spectrum including multiplicities. 

Example 6. Let K n = (V, E, uj) be the complete graph on n vertices with unit edge 
weights. That is, V = {vi, . . . ,v n } and 

w(eij) = < 1 J . for 1 <i,j <n. 
I otherwise 

(See figure^ left for n — 6. J 

For any fixed £ V note that V — {«&} 6 st(K n ). This follows from the fact 
that condition (i) of deHnition \3.2\ is trivially satisfied by the vertex set V ~ {vk}- 
Condition (ii) follows from the fact that ui(ekk) = ^ A. Therefore, the reduced 
graph TZy—{ Vk }(K n ) = (V — {«&},£ is well defined having edges e-ij with edges 
weights 

( \ Jl + l/A . 
(en = < tor 1 < i, 7 < n. 

I otherwise 

(See figure^ right for n = 6.) 

To compute the spectrum o~(K n ) note the following. As the matrix M{K n ) + / is 
the nxn matrix of ones then a(M (K n ) + 1) = {n, 0, . . . , 0} where has multiplicity 
n — 1 . Therefore, the graph K n has spectrum o~(K n ) = {n — 1, — f , . . . , — f } where 
— 1 has multiplicity n — 1. In particular, ^ o~(K n ). 

Via equation ([?}) the sets a(K n \{ Vk y) = {0} and <j~ 1 (K n \{ Vk )) = since the 
weight uj(ekk) — 0. As o~(K n \{ Vk y) n o~(K n ) — then theorem \3.5\ (or corollary^) 
implies a(K n ) — o~(7ty_t Vk y(K n )). That is, no eigenvalues are lost in reducing the 
graph K n to the graph TZ v _^ Vk y(K n ). 

We note that both <j(G\g) and <7 _1 (G|g) are easily calculated via equation (|7J|. 
Therefore, theorem offers a quick way of computing the eigenvalues of a reduced 
graph if the spectrum of the original unreduced graph is known. 

More importantly, however, theorem 13 . 5 1 tells us the eigenvalues we will gain and 
possibly lose if we reduce a graph. In section 6 this will allow us to modify the 
dynamical network (J 7 , X) while maintaining the network's stability. 
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k 6 n s (K 6 ) 

a(K 6 ) = <j(1Is(K 6 )) 

Figure 4. Reduction of K 6 over S = {vi, . . . ,v 5 }. The edges of 
K e and TZs{K 6 ) have weight 1 and 1 + 1/ A respectively. 



3.3. Sequential Reductions. In section 3.2 we observed that any reduction lZs{G) 
of a graph G G G is again a graph in G. It is therefore natural to consider sequential 
reductions of a graph G G G. This requires that we hrst extend our notation to 
sequences of isospectral reductions. 

For G = (V,E,lj) suppose S m C S m -i C • • • C Si C V such that Si G st(G), 
K 1 (G)=K Sl (G) and 

S l+1 G st(Ki(G)) where ft s<+1 (7^(G)) = 1 < i < m- 1. 

If this is the case we say Si, ... , S m induces a sequence of reductions on G with 
final vertex set S m . By way of notation we write lZ m (G) = 1Z(G; Si, . . . , S m ) where 
1Z(G; Si, . . . ,S m ) denotes the graph G reduced over the vertex set Si then S*2 and 
so on until G is reduced over the final vertex set S m . 

Theorem 3.6. For G = (V,E,u) suppose S G st(G). If there are vertex sets 
S C S m -i C • • • C Si C V then Si, ... , SVn-i, S induces a sequence of reductions 
on G. Moreover, 1Z(G; Si, ... , SVn-i, S) — lZs(G). 

Let S be a structural set of the graph G = (V, E, uj). In light of theorem 13.61 any 
sequence Si, ... , S m ~i, S where 

(8) S C S m -i C-CS 1 C7 

induces a sequence of reductions on G. To see how many possible sequences satisfy 
© we note the following. 

Each such sequence satisfying © corresponds to a partition of the set V — S. 
To see this suppose V = So and S — S m . Then the sequence Si, . . . , S m —i, S m 
corresponds to the partition {Pi, . , . , P m } of V — S where Pi — Si- 1 — Si for 
1 < i < m. 

Note that the number of ways to partition a set is exponential in its number of 
elements. Hence, there are often a large number of ways to sequentially reduce a 
graph to one of its structural sets. However, despite the potentially large number of 
such sequential reductions the following corollary of theorem [3lj] offers the following 
uniqueness result. 

Corollary 4. Suppose G — {V, E,uj) has structural set S. If S C S m -i C • • • C 
Si C V and S C T„_x C ■ • • C Ti C V then it follows from lemma [J7E\ that 

1Z(G; Si,..., S m -i, S) = TZs(G) = TZ(G; T\, . . . , T n -i, S). 

That is, the graph lZs(G) resulting from a sequence of reductions depends only 
on the final vertex set S if S G st{G). 
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ft(G;Si) K(G;Si,S 2 ) = n S (G) 

Figure 5. Sequential reduction of the graph G from figure [H 

Example 7. Let G again be the graph considered in example [^} If we let both 

51 = {vi, V2, V3, V4,} and S2 = {viyV^} then the following holds. 

The vertex set S2 is a structural set of G. In particular S2 = S where S is 
the structural set considered in example^ Moreover, S2 C Si. Hence, lemma 
implies the sequence Si,S 2 induces a sequence of reductions on G. Additionally, as 

52 = S then 1Z(G; Si, S2) = TZs(G) (see figures^ and\?$. 

We note that as S has four elements then there are fourteen ways to reduce G 
to the graph lZs{G). This follows from the fact that there are fourteen ways to 
partition a set with four elements. However, each of these reductions will result in 
the same graph. 

Let G = (V, E, oj) and suppose S £ st(G). It is a natural question to ask whether 
the result(s) of lemma [3T6l holds for the vertex set S. Specifically, (i) do any sets 
satisfying S C S m -i C ■ • • C Si C V induce a sequence of reductions on G; and (ii) 
can G be reduced to a graph with vertex set S via some sequence of reductions? If 
so, is this reduction unique? 

The answer to (i) is negative since Si need not be a structural set of G. To 
answer question (ii) we note the following. If the weight uj(eu) ^ A for some Uj G V 
then the vertex set S = V — is a structural set of G. This follows from the fact 
that S = {vi}. Hence, the graph G\§ is the graph restricted to the single vertex m. 
In particular, this implies that any cycle of G\§ is a loop. 

Therefore, any graph GeG can be reduced over the structural set S — V — {vi} 
if it is known that u>(ea) ^ A. Another way to state this is that it is possible 
to remove the vertex Vi from G via an isospectral reduction if tu(eu) ^ A without 
knowing anything about the graph structure of G. This has the following important 
implication. 

Suppose it is known that no loop of G or any loop of any sequential reduction 
of G has weight A. If this is the case then it is possible to remove any sequence of 
single vertices from G via a sequence of isospectral reductions. Therefore, G can 
be sequentially reduced to a graph on any subset of its vertex set. This idea is the 
motivation behind the following. 

For any polynomial p £ C[A] let deg(p) denote the degree of p. If w — p/q G W[A] 
where p,q G C[A] let 

ir(w) = deg(p) - deg(q). 
Again it is required to show that the function ir : W[A] — > Z is well defined. 
To do so suppose p/q and r/s are equivalent in W[A] i.e. ps = rq. We then have 
that deg(ps) = deg(rq). As deg{ps) — deg(p) + deg(s) and similarly deg(rq) = 
deg(r) + deg(q) then 

deg(p) - deg(q) = deg(r) - deg{s). 
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Figure 6. Distinct sequences of isospectral reductions with the 
same outcome. 



That is, ir(p/q) = n(r/s) implying the function n is well defined. 
Let W^A] be the subset of W[A] given by 

W W [A] = {we W[A] : 7t(w) < 0}. 

That is, W„-[A] is the set of rational functions in which the degree of the numerator 
is less than or equal to the degree of the denominator. Let be the graphs in G 
with edge weights in the set W T [A]. 

Lemma 3.7. If G G G w and S G st(G) then lZs{G) G G^. In particular, no loop 
of G and no loop of any reduction of G can have weight A. 

By the reasoning above, if G G G T then G can be (sequentially) reduced to a 
graph on any subset of its vertex set. This result is stated in the following theorem. 

Theorem 3.8. (Existence of Isospectral Reductions Over any Vertex Set) 

Let G = (V, E,uj) be graph in G OT and suppose V is a nonempty subset of V . Then 
there exist sets V G S m —i C • • • C Si C V such that Si, . . . ,S m -\, V induces a 
sequence of reductions on G. 

For G G it is therefore possible to reduce a graph G G to a graph on 
any (nonempty) subset of vertex set via some sequence of isospectral reductions. 
Moreover, such sequences have the following uniqueness property. 

Theorem 3.9. (Uniqueness of Isospectral Reductions Over any Vertex 
Set) Let G = (V,E,u>) be graph in G^ and suppose V is a nonempty subset of V . 
If Si, . . . , S m -i, V and Ti, . . . , T„_i, V both induce a sequence of reductions on G 
then K(G; S u . . . , 5 m _i, V) = K(G; T u . . . , T„_ 1; V). 

The results of theorem 13.81 and theorem 13.91 allows us to give the following defi- 
nition. 

Definition 3.10. Let G = (V, E,u>) be graph in G ff . If V C V is nonempty define 

K v [G]=1l(G;S l ,...,S m -i,V) 

where Si, ... , S m -i, V is any sequence that induces a sequence of reductions on G 
with final vertex set V. 
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Figure 7. G and H are equivalent under the relation induced by 
the rule r given in example |9] 

The graph 1Z\> [G] is well defined as a result of theorems l3.8l and l3.9l The notation 
IZv [G] given in definition 13.101 is intended to emphasize the fact that V need not 
be a structural set of G. 

Remark 7. Note that tt(c) = for any c e C. Hence, if M(G) e C" x " then 
G € Gtt . Therefore, any graph with complex weights can be uniquely reduced to a 
graph on any nonempty subset of its vertex set. This is of particular importance 
for the estimation of spectra of matrices with complex entries in |10) . 

Example 8. Let G — (V, E, lo) be the graph shown in figure® Our goal is to reduce 
G over its vertex set {v±, v^} C V ■ Note that as G S G w theorem ] 3.8\ guarantees that 
there is at least one sequence of reductions that reduces G to the graph 1Zi Vl!V4 \[G]. 

In fact there are exactly two. This follows from the fact that {vi,v^} ^ st(G). 
Hence, G cannot be reduced over {v\,Vi} with a single reduction. However, any 
(nontrivial) reduction of G removes at least one vertex from G. 

Therefore, the two possible ways of reducing G to the vertex set {^1,^4} are 

(9) TZ{ VljVi y[G]=TZ(G;{v 1 ,v 2 ,v i },{v 1 ,v i }); and 

(10) Tl{ VliVi }[G] = TZ(G; {vi,v 3) u 4 }, {v 1 ,v 4 })- 

Both of the reductions given in 0) and U0\) are shown in figure [H The dashed 
arrows labeled TZt in this figure represent the reduction of a graph over some struc- 
tural set H . This notation is meant to emphasize that this diagram commutes. That 
is, 

7^-{v 1 ,V4}(T^-{v 1 ,V2,V4,}{G)^) 1Z{ Vi Vi y(lZ{ Vl V3 Vi }(G)j 

as guaranteed by theorem \3.9\ . 

3.4. Equivalence Relations. Theorem 13.81 and theorem 13.91 assert that a graph 
G € G w has a unique reduction to any (nonempty) subset of its vertex set via some 
sequence of isospectral reductions. In this section this property will allow us to 
define various equivalence relations on the graphs in G n — 0. 

Two weighted digraphs G\ — (Vi, -Ei, wi), and G2 — (V2, -E2J 0^2) are isomorphic 
if there is a bijection b : Vi — > V2 such that there is an edge in G\ from v,i to Vj if 
and only if there is an edge e,j between b(vi) and b(Vj) in Gi with u>i(fiij) = wi(ey). 
If the map b exists it is called an isomorphism and we write G\ ~ G2 ■ 

An isomorphism is essentially a relabeling of the vertices of a graph. Therefore, 
if two graphs are isomorphic then their spectra are identical. 
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G H K 

Figure 8. G = U S [H] and K T [H] = K for S = {v l ,v 2 } and 
T = {v3,V4} but the graphs G and K do not have isomorphic 
reductions. 



Theorem 3.11. (Spectral Equivalence) Suppose for any graph G = (V,E,ui) 
in &„■ — that r is a rule that selects a unique nonempty subset r(G) C V . Then 
t induces an equivalence relation ~ on the set G„ — where G ~ H if the graph 

n T{G) [G]~n r{H) [H]. 

Example 9. Let G = (V,E,u)). The out degree of a vertex Vi G V is the number 
of outgoing edges incident to vi i.e. the number \{j : w(eij) ^ 0}|. If G G G^ Zei 
r(G) QV be the set of vertices of maximal out degree. 

Observe that for each graph G G G^ — the set t(G) both exists and is unique. 
Thus the relation of having an isomorphic reduction with respect to this rule induces 
an equivalence relation on — 0. 

In figure^ the graphs G and H have the vertex set t(G) — {fi,^} = t(H) 
of maximal out degree. As shown in the figure, the graph 1Z T ra\[G] ~ TZ t (h)[H]. 
Hence, G ~ H under the relation ~ induced by the rule r. 

Note that the relation of simply having isomorphic reductions is not transitive. 
That is, if 1Zs[G] ~ TZt[H] and Hu[H] ~ lZv[K] it is not necessarily the case that 
there are sets X and Y, subsets of the vertex sets of G and K respectively, such 
that TZ X [G] ~ n Y [K}. 

As an example, in figure |8] both TZs[G] ~ TZs[H] and TZt[H] ~ 7^t[^"] where 
5 = {^1,^2} and T = {v^, V4}. However, one can quickly check that for no subsets 
ICS and Y C T are lZx[G] ~ Overcoming this intransitivity requires 

some rule r that selects a unique set of vertices from each graph in G (see theorem 
I3.11|) . If r does not, e.g. r :=remove any vertex, then r does not induce an 
equivalence relation on G — 0. 

Importantly, choosing a rule that does select a unique vertex set allows one to 
study the graphs in G^ — modulo some particular graph feature. For instance, 
in example [9] this graph feature (or vertex set) is the vertices that do not have 
maximal out degree. 

4. ISOSPECTRAL TRANSFORMATIONS OVER FIXED WEIGHT SETS 

The isospectral graph reductions of the previous section modify not only the 
graph structure but also the weight set of a graph. That is, if lZs(G) — (S, £, fi) is 
any reduction of G = (V,E,lu) then typically w(E) ^ £«(£), i.e. 

{w(ejj) : etj G E} ^ {/i(ey) : G £}. 

This may lead one to assume that our procedure simply shifts the complexity of 
the graph's structure to its set of edge weights. However, this is not the case. 
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In this section we introduce different procedures of reducing and expanding a 
graph that restrict the weights of the transformed graph's edges to particular sub- 
sets of W[A]. As before, the procedure preserves the spectrum of the graph up to 
a known set. Such transformations are of particular importance in section 6 where 
dynamical network expansions are discussed. 

4.1. Branch Expansions. The idea behind an isospectral graph transformation 
that preserves a graph's edge weights is simple enough. If two graphs G,H E G 
have the same branch structure (including weights) then they should have similar 
spectra. If this is the case, then we say that G is an isospectral transformation of 
H over its weight set and vice-versa. 

To make this precise suppose G = (V,E,uj) and S E st(G). If the branch 
(3 = v il , . . . , v im G Bs(G) let Qg(P) be the ordered sequence 

CI G (I3) = oj(e ill2 ), . . .,u{e ij _ ui] ),uj{e ljlj ),uj{e lj ^ j+1 ), . . . , w(e im _ ljim ). 

for to > 1 and u^e^^) if to = 1. We call Qg(/3) the weight sequence of the branch /3. 
Moreover, we let ila(P)jj+i = u(e ijtij+1 ) and Q,G{P)jj = u(e iji:j ) f° r 1 < 3 < m-1. 

Let G, H E G. Suppose S — {v±, . . . , v m } is a structural set of both G and H. 
The branch set Bij{G; S) is isomorphic to Bij(H; S) if there is a bijection 

b:B ij (G;S)^B ij (H;S) 

such that Og(/3) = Qij(6(/3)) for each /3 e Bij(G; S). If such a map exists we write 
B ij (G;S)~B ij (H;S). If 

Bij{G; S) ~ Bij(H; S) for each 1 < ij < m 

we say Bs(G) is isomorphic to Bs{H) and write Bs(G) ~ Bs{H). 

Example 10. Suppose H and G are the graphs in figure^ with unit edge weights. 
The vertex set S = {vi,Vq} can be seen to be a structural set of both H and G. 
Moreover, 

Bn(H;S) = {vx,v 2 ,vi}, Bn(G;S) = {vi,u 5 ,vi}', 
B 13 (H;S) = {v!,v 2 ,v 3 }, B 13 {G;S) = 1*2,^3}; 
B 3 i(H;S) = {w 3 ,W4,wi}, B 31 (G;S) = {wi,u 4 ,u 3 }; 
B 33 (H;S) = {v 3 ,v 3 ;v 3 ,V4,v 3 }, B 33 (G; S) = {v 3 , v 3 ; v 3 , v u , v 3 }. 

Note that the branch ft — {vi, U5, vi} in B±i(G; S) has the weight sequence given 
by J7g(/3) = 1,1,1- Similarly, the branch 7 — {vi,V2,vi} in Bn(H;S) has the 
weight sequence f2#(7) = 1,1,1. Hence, B\\{H;S) ~ Bn(G]S). Continuing in 
this manner, one can check that each Bij(H;S) ~ Bij(G;S) for each i,j G {1,3}. 
Therefore, B S {H) ~ B S {G). 

Remark 8. We note that if Bs(H) — Bs{G) the graphs G and H need not be 
isomorphic (see example I10[) . 

For G G G and S G st(G) suppose a — Vx,...,v m and j3 — Ui,...,u n are 
branches in Bs(G). These branches are said to be independent if 

{u 2 , ■ • ■ ,v rn -i} U {u 2 , ■ . . ,m„-i} 7^ 0. 
That is, a and /? are independent if they share no interior vertices. 
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Definition 4.1. Let G,H eG and S G st(G), st{H). Suppose 

(i) B S (G) ~ Bt(H); 

(ii) the branches of Bs{H) are independent; 

(hi) each vertex of G and H belongs to a branch of Bs{G) and Bs(H) respectively. 
Then we call H a branch expansion of G with respect to S. 

Proposition 1. Let GsG and S G st(G). If H and K are branch expansions of 
G with respect to S then H ~ K . 

Proof. Suppose H — (Vi,£?i,a;i) and K = (V2, E%, W2) are branch expansions of G 
with respect to S 1 . Then by assumption (i) of definition 14.11 it follows that there is 
a bijection 

b : B S (H) -> B S (G) 
such that Q H (/3) = ft x (&(/?)) for all /3 G £<?(#). 

For the branch f3 = , . . . , Vi m let = Wj . for 1 < j < m. We call j the 
index of Vi j . By assumption each vertex v G V\ belongs to a branch of Bs{H) and 
the branches of Bs{H) are pairwise independent. Hence, either v G S or v — f3(i) 
for exactly one (3 G Bs(H) and index i. Let £? : V\ — > V% be the map 



B(v) 



v if v G S 1 

(6(/3))(i) if w £ 5 and v = (3(i) 



By assumption each v G V2 is similarly either in 5 1 or equal to /3(i) again for 
exactly one /3 G Bs(K) and index i. Therefore, B has the inverse 



B~\v) 



v if v G S 

(6 -1 ()8))(i) if w ^ 5 and « = /8(i) 

implying the map B : Vi — >■ V2 is a bijection. 

Let ey G -Ei. If t^, Vj G S 1 then Wj, Vj G Bs{H) implying 

wi(eij) = Vt H (vi,Vj) = VL K (v B{€jl v B{j) ) = u>2(e B {i)B{j))- 

Suppose Vi, Vj ^ S. If Vi ^ Vj then there are branches a\ — v\, . . . , Vi, . . . , v s and 
Q2 = v>i, . . . , Vj, . . . , Vt in Bs(H) containing Vi and Vj respectively. Therefore, the 
sequence V\, . . . , Vi, Vj, . . . , Uf G Bs{H). As the branches of Bs{H) are pairwise 
independent then ot\ = a.2 and j = i + 1. Hence, 

Wl(eij) = QH(ctl)i,i+l = £lK(K a l))i,i+l = ^2{dB(i)B( 3 ))- 

If Vi = Vj then there is a unique /3 G Bs{H) and index k such that = /3(h). 
Hence, 

wi(eij-) = SlH{P)kk = &K{b{fi))kk = uj 2 {eB(i)B(]))- 
Supposing Vi G S, Vj ^ S or Wj ^ 5, Vj G 5 a similar argument implies that 
wi(ey) = w 2 (e B (i)B(j))- Hence, H ~ K. □ 

It follows that a branch expansion of G with respect to S is unique up to a 
labeling of vertices. Therefore, any two expansions of G with respect to 5* are 
isomorphic. In what follows we let Xs{G) be a representative of the class of branch 
expansions. By slightly abusing our terminology we call any representative Xg{G) 
the branch expansion of G with respect to S. 

The principle idea behind a branch expansion is the following. If G G G and 
S G st(G) then the set of branches Bs{G) is uniquely defined. However, there are 



ISOSPECTRAL TRANSFORMS, SPECTRAL EQUIVALENCE, AND GLOBAL STABILITY 23 




Figure 9. The graph G is a branch expansion of the graph H 
with respect to the structural set S — {vi,vs}. 



typically many other graphs H with the same branch structure as G, i.e. S £ st(H) 
such that B S {H) ~ B S (G). 

A branch expansion of G over S is then a graph H = Xs(G) with identical 
branch structure but with the following restriction: The branches of Bs(H) are 
pairwise independent and every vertex of H belongs to a branch in Bs(H). That 
is, any vertex of S in H is part of exactly one branch in Bs(H). 

Hence, given a graph G and structural set S we can algorithmically construct the 
expansion Xs(G) as follows. Start with the vertices S. If /3 € Bij(G; S) then both 
Vi,Vj £ S. Construct a path (or cycle) from Vi to Vj with weight sequence f2(/3) 
with new interior vertices. By new we mean vertices that do not already appear on 
the graph we are constructing. Repeat this for each j3 € Bij (G; S) . The resulting 
graph is the branch expansion Xs{G). A branch expansion is our first example of 
an isospectral graph transformation that preserves the weight set of a graph. 

Theorem 4.2. Let G = (V,E,lo) with structural set S. Then the graph G and its 
branch expansion X$(G) have the same set of edge weights. Moreover, 

det (M(X S (G)) - XI) = det (M(G) - XI) J] (cu(e u ) - A)™ ; 1 
where rii is the number of branches in Bs(G) containing Vi. 

Example 11. Consider the graph G — (V, E, lu) and H — (V, £ , /i) in figure^ with 
unit edge weights. As demonstrated in example \10\ if S = {wi,^} then the branch 
set Bs{G) ~ Bs(H). Moreover, it can be seen from figure^ the five branches of 
Bs(G) share interior vertices. That is, the branches of Bs(G) are pairwise inde- 
pendent. Lastly, each vertex of G belongs to at least one branch ofBs(G). 

Therefore, the graph G is a branch expansion of H with respect to S, i.e. G = 
Xs{H). Importantly, note that the edge weights of H and its expansion G are 
identical, i.e. the sets uj{E) and n(£) are both {0,1}. 

Observe that the vertex V2 of H is an interior vertex of the branches V\ , v% , V% ; 
^i , V2 , W3 £ Bs {H) . Similarly, the vertex V4 of H is an interior vertex of the 
branches W3,W4,W3; vs,V4,vi € Bs(H). As uj{e22) = 1 and w(e44) = theorem 
\4l\ implies that a(X s (H)) = a(H) U {1, 0}. 

We note that G — X$ (H) in figure is isomorphic to the graph G in figure [H 
Hence, a{G) = {2,-1,1,1,0,0} implying a{H) = {2,-1,1,0}. 

4.2. Isospectral Graph Transformations Over Modified Weight Sets. The 

branch expansions of the previous section are a method of graph transformation 
that separates the various branches of a graph. In this section we consider the 
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reversal of this process. Specifically, we introduce a method of transforming a 
graph that incorporates a new technique of merging branches. 

This type of isospectral transformation will have the additional property that 
it keeps the graph's weights in a fixed subset U C W[A]. The particular type of 
subsets for which this will hold are semirings of W. 

The set U is a semiring of W[A] if it has the following properties. Both 0,1 £ U 
and for any u\,U2 £ U both the product u\U2 and the sum u± + U2 are in U. 
Moreover, as U C W[A] we are implicitly assuming that if u € U then any other 
representation of u is also in U. Examples of such semirings include C[A], K, and 
Z+ = {0,1,2,...}. 

Definition 4.3. For G = (V, E,ui) let S £ st(G). The set S is a complete structural 
set of G if 

(i) each cycle of G, including loops, contains a vertex in S; and 

(ii) ui(eu) ^ A for each Vi £ S. 

The difference between a structural set and a complete structural set of a graph 
G is the following. Recall that if S is a complete structural set of G then every 
cycle of G contains a vertex in S. If S is simply a structural set of G then loops of 
G need not contain a vertex of S. 

For G £ G let sto(G) denote the set of all complete structural sets of G. Addi- 
tionally, let (7o (G) be the nonzero elements of o~(G) including multiplicities. Hence, 
we refer to o~q(G) as the nonzero spectrum of G. Moreover, we denote by p(G) the 
spectral radius of G. That is, 

piG) = max |4 



The following is a corollary of theorem [ 
Corollary 5. Let G = (V,E,lj) and suppose S £ sto(G). Then 

a (U s (G)) =* (G). 
In particular, p(lZs{G)) = p(G). 

Proof. Suppose G = (V, E,ui) and S € sto(G). As no vertex in S can have a loop 
then u!(eu) = for each Vi € S. Equation ([7]) then implies that 

det(M(Gy - XI) = J[ A. 

Hence, a- 1 (G\ s ) = and a(G\ s ) = {0, ... ,0} in which has multiplicity The 
corollary then follows from theorem 13.51 □ 



To simplify the discussion in what follows, suppose G 6 G and S E st(G). If 
j3 = ?;,!,... , Vi m £ Bs(G) we will say ei j .i j+1 is the jth edge belonging to the branch 
(3 for each 1 < i < m — 1. 

Lemma 4.4. Let G = (V,E,lj) and S £ st (G). Then X S {G) = (V,£,p) has the 
following properties. 

(i) If eij £ £ then eij belongs to exactly one branch of Bs{Xs(G)). 

(ii) If ft = v\, . . . , v m is a branch in Bs(Xs(G)) then 

tt Xs (G){0) = M(ei2), • ■ ■ , v(ek,k-i), Q, p(ek,k+i), M e m-i,™)- 
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Proof. As each vertex of Xs(G) belongs to at least one branch of Bs(Xs{G)) then 
the same holds any £ £. On the other hand, suppose dj belongs to both 
/?i,/?2 £ Bs(Xs(Gj). Then neither Vi or Vj can be interior vertices of (3\ or /3 2 as 
these branches, if distinct, are independent. Hence, /3± = = Vi,Vj or ey belongs 
to at most one branch of Bs(Xs(G)). This verifies property (i) 

Since S is a complete structural set of G then oj(eu) = for each Vi £ V — S. 
This implies that, for any /3 = Vi, ...,v m £ Bs(G) the weight sequence flc(P) has 
the form 

Q.g(/3) = w(ei 2 ), • ■ • ,u(e k ,k-i),0,u(e k , k+ i), . . . ,w(e TO _i )m ). 
As B S (G) ~ 6 s (,*s(G)) property (ii) holds. □ 

The ^-Construction: Branch Reweighting 

Given a branch expansion X$ (G) we construct a new graph 3^s (G) by reweighting 
the branches Bs{Xs{G)) of Xs(G). The idea behind this construction is to reweight 
the branches Bs{Xs(G)) in such a way that it preserves their branch products. 

Suppose the expansion Xs(G) = (V,£,[i) where S £ sto(G). Let the graph 
ys(G) = (y,£,v). That is, ys(G) has the same vertex and edge set as Xg(G) 
but possibly different edge weights. We note this implies that S is a complete 
structural set of ys(G) and moreover that the branch set Bs(ys(G)) is identical 
to B S (X S (G)). 

For the branch (3 — v\, . . . ,v m £ Bs(ys(G)) define the weight sequence of /3 to 

be 

7YI — 1 

(11) ^ s (G)(/3)= II M(efc,fc+i),0,..., 1,0,1,.- -,0,1 

fe=i 

if to > 1. If to = 1 let fi;y s ( G )(/3) = ^i(en). As Xs(G) and ys(G) have the same 
vertex and edge set then lemma implies that each edge of 3^s(G) belongs to 
exactly one branch of Bs(ys(G))- Therefore, equation (fTTj) completely specifies 
the edge weights of the graph 3%(G). 

In particular, the edge £ £ in ys(G) has weight v(e.ij) = 1 unless is the 
first edge of a branch in Bs(ys(G))- If happens to be the first edge of the branch 
P G Bs(ys(G)) then its weight is the product of the nonzero entries of Qx s (G)(l3)- 

Considering part (ii) of lemma l4.4[ 3^s(G) is effectively the graph X$(G) in which 
the first edge weight of each branch is the product of that branches original weights 
in Xs(G). Every other edge of ys{G) is given weight 1. fFigurellOl gives an example 
of this branch reweighting.) 

The Z-Construction: Branch Merging 

From the branch reweighting ys(G) we construct the graph Zs(G). The major 
idea behind this construction is that the branches Bs(ys (G)) of the graph ys(G) 
can be merged together in a way that maintains the weight set of the graph as well 
as its nonzero spectrum. 

For G = (V, E, uj) suppose S — {v%, . . . , v m } is a structural set of G. Let 

Bj(G;S)= |J B^G-S). 

l<i<m 

That is, Bj(G; S) are the branches in Bs{G) terminating at the vertex vj. For any 
/3 = U\, . . . , u k in Bs{G) let \(3\ = k, i.e. the number of vertices in /3. Moreover, let 
{(3}int denote the set of interior vertices of j3. 
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Suppose the graph ^s(G) = (V, £ , v) and S — {vx, . . . ,v m } 6 sto(G). For each 
Vj £ S select a branch ft e Bj(ys(G)) with the property that |^'| > \ft for all 
branches /3 G Ss(^s(G)). We set ft = if Bj(^ s (G)) = 0. 

Denoting B = B 3 (y S (Cf)) ~ {ft, ■ ■ • , P m ] let 

(12) W = (J {/3} mt . 

That is, U is the set of interior vertices of the branches (3 ^ {ft, . . . , f3 m } in 
B s (y s (G)). Let the graph 

Z' s (G)=y s (G)\ V -u. 

Recall that the branches of Bs{ys{G)) are pairwise independent. Therefore, the 
graph Z' S {G) = (V — U,£' ,v') where e € £' if and only if e belongs to some ft . 
Furthermore, the edge weights of Z' S {G) are given by the restriction v' = v\s>. 

If e is an edge from the vertex a to the vertex b we will denote this by e = (a, b). 
Suppose the branch ft = u{ , . . . , v k . For each 

(3eB l3 (y s (G);S)-ft 

we add an edge (Vi,v J k _i^, +2 ) to the graph Z' S (G). The edge {vi,v 3 k _^^ +2 ) is given 
the weight of the first edge belonging to (3 in ys(G). If this is done over all 
1 < i,j < m we call the resulting graph Zg(G). 

It is important to note that the graph Zg{G) may have parallel edges. By parallel 
edges we mean that there may be multiple edges in the edges set of Zg(G) of the 
form (a,b). In particular, if two branches ft, ft 6 Bij(ys{G); S) have the same 
length, i.e. \ft \ = \ft\ = £, then there are (at least) two edges in Zg{G) of the form 

( u i> ^-|^|+2)' 

Suppose Zg(G) has parallel edges ei,...,ejv of the form {v^Vj) with weights 
wi, . . . ,wjy. We replace the edges ex, ■ ■ ■ , &n in Zg(G) with the single edge eij 
having weight wi + ■ ■ ■ + . If this is done for each set of parallel edges in Zg(G) 
we denote the resulting graph by Zs(G). 

Note that our construction of Z${G) depends on the initial choice of each ft. 
We therefore write Z S (G) = Z S (G; ft, . . . , j3 m ). 

Theorem 4.5. Let G — (V,E,lu) and S G s£o(G). Suppose the edge weights of 
G are in the semiring U C W[A]. Then Zs(G) has edge weights in U. Moreover, 
a Q (Z s (G)) = a (G) " ' 

Example 12. Let G = (V, E, w) be the graph shown in fiaure [W\ (far left). Note that 
the vertex set S — {vi, ^2,^3} is a complete structural set of G since G has no cycles 
(including loops). The branches of G with respect to S are then ft — vi,V3,Ve, 
ft = v\, V4, vq, ft = vi, V4, v-j, and ft = Vi,V2,Vs,V7. Observe that only ft and ft 
share an interior vertex. 

The graph Xg{G), shown in haure XlXh (middle left), is then the graph in which the 
branches ft and ft have been replaced with the independent branches ft = vi,ui, vq 
and ft = v\,U2,vj respectively. Thus, the expansion Xs{G) has the vertex set 
V = {vi,v 2 ,v 3 ,v 5 ,v e ,v 7 ,u 1 ,u 2 }. 

Note that the products of the weights along the branches ft, ft, ft, ft are 
5, 12, 14, 96 respectively. Hence, these are the first edge weights of each of the 
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Vi Vi Vi Vi 



G X S (G) y s (G) Z S {G) 

Figure 10. The isospectral transformation Zs(G) of G over S — {vi, 1)%, V3}. 

branches Pi, p 2) Pa, Pa in ys(G) respectively (see figure \T0\ middle right). Each 
other edge weight ofys{G) is 1. 

To construct Z' S (G) note that P 1 = 0, P 7 = Pa, and P 6 is either Pi or p 2 . Here 
we make the arbitrary choice of letting p e — Pi- As the sets {ftjint = {ui} and 
{Pz]int = {u 2 } then Z' S {G) has vertex set V - {ui,u 2 } = {vi,v 2 ,v 3 ,v 5 ,v 6 ,vr}. 
Since the branch set 

B s (Z> s (G))-{p e ,p 7 } = {p 2 J 3 } 

then constructing Zg(G) amounts to adding two edges to Z' S (G); one for P2 and 
one for P3 . 

Note that p 2 E B 16 (y s (G);S) with first edge weight 12 and \p~ 2 \ = 3. As \P e \ = 3 
then to the graph Z' S (G) we add a parallel edge from vi to U3 with weight 12. For 
the branch P% G Bn(ys(G); S) note that it has first edge weight 14 and \p 2 \ = 3. 
Since \p 7 \ = 4 then to the graph Z' S (G) we add an edge from vi to v$ with weight 
14. The result is the graph Zg{G). 

Note that the two parallel edges from v\ to V3 in Zg(G) have weights 5 and 
12 respectively. Hence, in Zs{G) = Zs(G] P\, Pa) shown in figure [771 (far right) 
the edge en has weight 17. Moreover, the edge weights of G are in the semiring of 
nonnegative integers Z + . Therefore, as guaranteed by theorem ^. 5\ the graph Zs(G) 
also has edge weights in Z + . One can also compute that ao(Zs{G) — ao(G). 

For a graph G = (V, E,ui) let \G\ = \V\, i.e. the number of vertices in V. 

Proposition 2. For G = {V,E,uj) let S — {v%, . . . ,v m } be a complete structural 
set of G. Then 

m 

\Z S (G)\ = |5| + ^(|^|-2) 
i=i 

where \P j \ = 2 if ft = 0. 

Proof. Let y s (G) = (V,U,v). Suppose a,P € B s {y S (G)) and a ^ p. Note that 
the branches of Bs (3^s (G) ) are pairwise independent and each vertex of V belongs to 
a branch of B s (ys(G)). Hence, {a} int n {P} ln t = and {P} int n S = 0. Therefore, 
V is the disjoint union 

V = SU( |J {Phut)- 
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Recall that the set IA 7 given by (|T2j) . is the set of interior vertices of the branches 
B s {y S {G)) -{Pi,...,Pm}. Thus, the vertex set 

m 

v-u = su(\J{ft} mt ). 

8=1 

where each {ft} lnt n {ft} ln t = and {ft} tnt n S = for i ^ j. Therefore, 

m 

|V -W| = |5| +£(|^'|-2) 
i=l 

since each branch /J- 7 ' has — 2 interior vertices. If ft — then it has no interior 
vertices which is compensated for by assuming \ft\ = 2. 

Note that \Z' S S(G)\ = \y S {G)\v-u\ = \V - U\. As constructing Zg(G) from 
Z' G (S) only involves the addition of edges and constructing Zs(G) from ZgG(S) 
the removal of parallel edges it follows that |Zs(G)| = |V — U\. This completes the 
proof. □ 

In example [T^] we consider the 2-construction of the graph G in figure [TU] over 
the complete structural set S — {v\,v&, v>j}. From proposition [2] it follows that 
\Z S {G)\ = \S\ + \ft\ + \ft\ + \ft\. As ft = 0, ft = v u v 3 ,v 6 , and ft = v 1 ,v 2 ,v 5 ,v 7 
then \Z S (G)\ = 6. 

Note that in this example \G\ > \Zs(G)\. Hence, Z${G) can be considered to be 
a reduction of G over the weight set Z + having the same nonzero spectrum. 

5. Proofs 

In this section we give proofs of the theorems and lemmas in sections 3 and 4. 
Let GeG. If S 6 st(G) then, as previously noted, the vertices of G can be ordered 
such that M(G) — XI has the block form 

' A B 



(13) M(G) -\I = 



C D 



where A corresponds to the vertices in S and is a triangular matrix with nonzero 
diagonal (see Frobenius normal form in [5]). The matrix A can be made triangular 
from the fact that G\§ contains no (nonloop) cycles. The nonzero diagonal follows 
from the assumption that uj(ea) =/= X for V{ G S. As a triangular matrix with 
nonzero diagonal has a nonzero determinant then the matrix A is invertable. 
Using the matrix identity 

(14) det ( C Dl = d0t ^ ' dct ^ ~ CA ~ lB ) 



it follows that 



6*{p-CA-iB) = ^ M{ ®-X) 



det(A) 

As D - CA^B = R- XI for some R e W[A]I 5 I X I 5 I then the claim is that the 
isospectral reduction TZg(G) is the graph with adjacency matrix R (see proof of 
lemma [5T3"|) . Furthermore, as A corresponds to the vertex set S then by equation 


det(A)= H (w(e«)-A). 

Vi<£S 

Our first task in this section is to make this precise. 
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Definition 5.1. Let M g W[A]" xn and IC{l,...,n} be nonempty. If the set 
X = {ii, . . . ,i m } where ij < ij+i then the matrix A given by Akt = Mi k i e for 
1 < k,£ < m < n is the principle submatrix of M indexed by X. 

In what follows, we will use the convention that if X = {ii, . . . , i. m } is an indexing 
set of a matrix then ij < ij+i- Moreover, if D the principle submatrix of the matrix 
M e W[A]" X ™ indexed by X then there exists a unique permutation matrix P such 
that 



where A is the principle submatrix of M indexed by X = {1, . . . , n} — X. 
Assuming A is invertable, then by (fT4|) 

det(M) = det(PMP- 1 ) = det(A) det(D - CA^B). 

If this is the case we call the matrix D — CA~ 1 B the reduction of M over X and 
write r(AI;X) = D — CA~ X B. In the case that X — {1, . . . ,n} we let the matrix 
r(M;X) — M. Moreover, if r(M;X) can be reduced over the index set J we write 
the reduction of r(M;I) over J by r(M;X, J), and so on. 

Lemma 5.2. Let M G W[A] nx ™ andX m C {1, . . . , n} be nonempty. If the principle 
submatrix of M indexed by X m is upper triangular with nonzero diagonal and the 
sets X m C Z m _i C ••■ C Xx C {l,...,n} then r(M;Xi, . . . ,I m _i,I m ) = r(M;X m ). 

Proof. Without loss in generality suppose X m C {1, . . . , n} such that the matrix M 
has block form 



where A is an upper triangular principle submatrix of M indexed by X m with 
nonzero diagonal. Then for X m C Ij C {1, . . . , n} there exists a unique permutation 
matrix P such that 



where Aq, A±, and D are the principle submatrices indexed by {1, . . . , n} — X±,X± — 
X mi and X m respectively. Moreover, both Aq and Ai are upper triangular matrices 
with nonzero diagonal as they are principle submatrices of A and are therefore 
invertable. 

By partitioning the matrix PMP^ 1 into the blocks 





PMP 



Aq B a B 2 
Co A l B l 
C 2 Ci D 




A 
Co 
Ci 



Bo B 2 
Ai B 1 
Ci D 




it then follows that 



(15) 




The claim is then that the matrix CqAq 1 Bo is an upper triangular matrix pos- 
sessing a zero diagonal. To see this let the sets {1, . . . , n} — X\ = . . . , i s }, 
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X\ — X m = {jx, . . . ,jt}, and (A ) = aij. As the inverse of an upper triangular 
matrix is upper triangular onj — for i > j. Hence, 

s p 
p=l g=l 

As 1 < q < p then each i q < i p . For each k < £ note that similarly jj~ < j(. 
Moreover, if < i p then A ip j k — since A is upper triangular. If j\ > i p then 
*g < ji implying Aj iiq = for the same reason. As jf. ^ i p , it then follows that 

(C A- x B ) lk = for all k < I 

verifying the claim. 

A\ — CqAq 1 Bq is therefore an upper triangular matrix with nonzero diagonal 
implying A\ — CoA ~ 1 Bo is invertible and r(M;Xi) can be reduced over X m . In 
particular, 

(16) r(M;Z 1 ,I m ) = D — C 2 A Q X B 2 - [C x - C 2 A Q 1 B )V(B 1 - C A^B 2 ) 

where V = (A 1 - C A^ 1 B )- 1 . 

To verify that r(M;Xi,X m ) = r(M;X m ) note that M has block form 



M 



A B 
C D 



and QAQ- 1 



Ao B 
Co Ai 



where Q is the principle submatrix of P indexed by X m . Therefore, 

(Ao - BoA^Co"- 1 a- 1 ^-! 



Q- 1 



1 


' Ao 


Bo ' 


-l 


r 1 




. C ° 


A i . 


Q = C 





-r-'CoA^ 1 



-Ao L B T- 



^o^o r 

Note the matrix (A - Bo^r^o) -1 = A^ 1 _ A Q 1 B Q TCoA 1 by the Woodbury 
matrix identity ( see [TB]) and is therefore well defined. From this 

n -1 

QB 



r{M-Xm) =D-CQ- 1 

=D-[C 2 Cx 





' B 2 ' 







CxTBx. 



Ao Bo 
Co A l 

(Ao - BoA^CoY 1 -A^BoT- 1 

-v-^CoA- 1 r- 1 

=D - C^A- 1 - A- 1 BoTC A- 1 )B 2 - CxTCqA^B^- 

C 2 A Q ~ 1 BqTCoAq 1 Bi 

From (Unj) h follows that r(M;X m ) = r(M;X 1 ,X m ) 

For X m CI 2 C X\ the same argument can be repeated to show 

r(r(M;X!);X m ) = r(r(M ; X x ); X 2 ; X m ) 

as r(M;Xi) has the same form as M, i.e. its upper left hand block is an upper tri- 
angular matrix with nonzero diagonal. Hence, r(M;X m ) = r(M;Xi,X 2 ;X m ). The 
lemma follows by further extending X m C X% C Ij C {1, . . . , n) to the nested 
sequence X rn C Xm-% C • • • C X± C {1, . . . , n} by repeated use of the above argu- 
ment. □ 

To establish that graph reductions are analogous to matrix reductions we prove 
the following lemma. 

Lemma 5.3. If S — {v m +\ 1 ■ ■ ■ ,v n } is a structural set of G G G then lZs(G) is 
the graph with adjacency matrix r(M(G) — XI; {m + 1, . . . , n}) + XI . 
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Proof. Let S — {t>i, . . . , v m }. From the discussion in the first paragraph of section 
5, if S is a structural set of G — (V, E, lu) then M(G) — XI has the block form 



M(G) - XI 



A B 
C D 



where A is the principle submatrix of M(G) — XI indexed by {1, . . . , m} and is an 
upper triangular matrix with nonzero diagonal. 

With this in mind, let Vk = {vk+i, ■ ■ ■ , v n \ and 1% = {fc+l, . . . , n} for 1 < k < m. 
If M(G)ij — u>ij for 1 < i,j i < n then, proceeding by induction, for k = 1 

r(M(G) - XI -li) = (Di - [«2i, . . . ,co nl f — r [wi3, . • .,«!„]) 

where D\ is the principle submatrix of M(G) — XI indexed by I\. Hence, 

(17) r(M{G)-XI-l 1 ) t] 



lo 1} + u>au>ij/(X - wn) i ^ j 
uiij + uiuuiij / (\ - u)n) -X i = j 



for 2 < i, j < n. 

Conversely, note that the set Bij (G; V\ ) consists of at most the two branches 
/3 + = Vi,Vi,Vj and f3~ = Vi,Vj i.e. the branches from Vi to Vj with and without 
interior vertex v\. As V u (/3~) — and 'P w (/3 + ) = UJnUij/ (X — oju) if the branches 
j3~ and (3 + are respectively in Bij(G\V\) it then follows from ^ and p7|) that 
r(M(G) - A/;X X ) = M(TZ Vl {G)) - XI. 

Suppose then that r(M(G) - XI]I k -\) = M('^v fc _ 1 (G)) - A/ for some k < m. 
Since the principle submatrix of M(G) — A/ indexed by Ik is upper triangular with 
nonzero diagonal and Ik C Zfc_i C {1, . . . , n} an application of lemma E21 implies 
that r(M(G) - XI:I k ) = r(M(G) - AJ;X fc _i,X fc ). Hence, 

r(M(G) - XI; I k ) = r(M(U Vk _ 1 {G)) - Wfe)- 
Letting = M(7^v fc _ 1 (G)) then by the argument above 



r(M(G) - XI:I k )ij = 



Mij + MikMkj I '(X - Mkk) i + j 

Mij + MikMkj I '(A - Mkk) - X i = j 



for k + 1 < i, j < n. As Mij — M(Rv k _ 1 (G))ij then the entries of M are given by 
Observe that 

E = E E 

/3eBy(G;V fe ) ^B-.(G;V k ) ,3G8+(G;V fe ) 

where B^(G;Vk) and S^(G;Vfc) are the branches in B,j(G; Vfc) that contain and 
do not contain the interior vertex Vk respectively. It then immediately follows that 
Bi 3 (G- V k ) = Bui?; V k -i) implying E^ B r.( G; v fc ) = ^.r 

On the other hand, any j3 G B^(G; Vk) can be written as /3 = Vi, . . . , Vk, ■ ■ ■ , Vj 
where f}\ = Vi,. . . £ S ifc (G; V fc _i) and /3 2 = Ufc,- • . ,«j G B k3 (G; V k ~i). Hence, 
equation ^ implies 

as) E w)= E g^Cfe) . 

X — LOkk 

0eB±(G;V k ) 0eB+(G;V k ) 
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Conversely, if /3i = v it . . . , Vk G B ik {G; V k -i) and f3 2 = v k ,..., vj G B k] (G; Vk-%) 
then /? = ^j, . . . , Vk, ■ ■ ■ , Vj G B^(G; V k ). This follows from the fact that /?i and 
/?2 share no interior vertices since otherwise G|y fc would contain a cycle. Note that 
G\y k cannot contain a cycle since V k G si(G). Therefore, 

(19) y, p«-ow»G8a)= E E 

/3ee+(G;V fc ) /3iei3 lfc (G;Vfc_i) /3 2 ee fc3 (G ; v fc _i) 

Moreover, Wfefe = M kk since B kk (G; V k -i) contains at most the cycle v k , v k given 
that V" fc _i G st(G). From (HHJ) and ([11]) it then follows that 

Mij + M lk M kj /(X - M kk ) = ]T P„C8) 

/3eBy(G;Vj«) 

implying r(M(G) - A/;Z fe ) = M(fty fc (G)) - XI. 

By induction we have that r(M(G) - A/;X m ) = M{TZ S (G)) - A/ by setting 
k = m. □ 

If S = {vi ± , . . . , Vi m } is a structural set of G € G then let {ii, . . . , i m } be the 
index set associated with S. Hence, if S € st(G) is indexed by I then S is indexed 
by X and there is a unique permutation matrix P such that 

A B 



P(M{G)-\I)P- 1 = 



C D 



where A is the principle submatrix of M{G) — XI indexed by X. Therefore, A = 
M(G\§) - XI and it follows from lemma [Ol and (fig]) that 

(20) det (M(G) - A/) = dct (M(G\§) - XI) det (M(^ S (G)) - A/) . 

Let det(M(G| s ) - XI) = p/q G W[A] where p,q G C[A]. Then theorem E3] follows 
by observing that a(G\g) and ct -1 (G|5) are the solutions to p = and = 
respectively. 

We now give proofs of the theorems on isospectral transformations. For theorem 
13.61 we give the following. 

Proof. Let S be a structural set of G = (V, E, lj). Following definition 13.21 the 
graph G\g contains no cycles, except possibly loops, and if G\§ has loops they 
do not have weight A. If 5 C T C V then the graph G\f is a subgraph of G\§. 
Therefore, T G st(G) as G\f similarly contains no cycles, except possibly loops, 
where loops do not have weight A. By the same argument, S G st(1ZT(G)) since 
H,t(G)\s is a subgraph of G\§. Hence, T, S induces a sequence of reductions on G. 

Suppose S m C S m -i C • • • C Si C V where S m is a structural set of G. It 
then follows that Si, ... , S m induce a sequence of reductions on G. Let It be the 
index set associated with Si for 1 < i < to. By lemma [5731 the graph TZ(G; Si) has 
adjacency matrix r(M(G) — XI, Xi) + XI. Hence, by another application of lemma 
15.31 the graph TZ(G; Si, S2) has adjacency matrix 

r(r(M(G) - AJ;^) + (A - X)I;X 2 ) + XI = r(M(G) - XI;X 1 ,X 2 ) + XI. 

Continuing in this manner, it follows that 1Z(G; Si, . . . , S m ) has adjacency matrix 
given by r(M(G) — XI;Xi, . . . ,X m ) + XI which by lemma [531 is equivalent to the 
matrix r(M(G) - XI;X m ) + XI. 

Given that r(M(G) — XI;X m ) + XI is the adjacency matrix of 7?-s m (G) then 
lemma IB~3l implies that TZs m (G) = TZ(G; Si, ... , S m ). □ 
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We now give a proof of lemma 13.71 

Proof. Let wi — px/qi and W2 = V2I12 be in W T [A] for f»i , P2 , 9i , 92 G C[A]. That 
is, 7r(w7i), k{w2) < 0. Then 

tt(wi + W2) = 7Tl I < max{7r(wij, ir{W2)} < 0, 

V 9l<72 ' 
/ P1P2 \ 

tt{wiW2) = 7r( = n(wi) + ir(u)2) < 0, and 

\ 9i 92 ' 



1 \ (1i 



-1 < 0. 



.X-WiJ v(9iA-pi)> 
It then follows from equation ([5]) and © that lZs(G) G G^ if G G G„-. □ 

In order to prove theorem 13.91 we require the following lemma that essentially 
implies that the order in which any two vertices are removed from a graph G G G n 
does not effect the final reduced graph. 

Lemma 5.4. Suppose G G G^ with vertex set V — {vi, . . . ,v n } for n > 2. // 
T = {v 1 ,v 2 } thenK(G;V - {v!},T) =K(G;V - {v 2 },T). 

Proof. Let G G Gtt and M(G%- = u l0 for 1 < i, j < n. Then M(G) - XI can be 
written in the block form 



cjn - A 

M(G) - XI = \ lj 2 i w 2 2 - A B 2 ) where A 




OJll — A LO12 
U)21 w 22 — A 



is the principle submatrix of M(G) — XI indexed by X — {1, 2}. 

Note the assumption G G G,r implies that both wu — A, CJ22 — A 7^ 0. Moreover, 
if det(A) = then CJ22 = ^\2^i\{^\\ — A) -1 + A. However, equation (jTTJ) implies in 
this case that 

M(TZ V - V1 (G)) 22 - A = L022 + - A = 0. 

Hence, M\R,y_ Vl (G)) 22 = A, which is not possible as it contradicts the conclusion 
of lemma 13.71 Therefore, det(A) ^ implying 

"" { D ~ L ~x ( ^ l! - - Srr ,:r ' ~ ^~x° 3m ~ ^~x Bj - 

is well defined. 

Letting Zj = I U {i} for i ^ j and i, j G I then repeated use of (fl4| implies 
= r(M(G) — XI;li,T). Moreover, as R12 = R21 then it follows that the graphs 
TZ(G; V — {vi}, T) = TZ(G; V — {^2}, T) since these graphs have the same adjacency 
matrices. □ 



To simplify the proof of theorem l3.9l we note the following. If Si, ... , S m induces 
a sequence of reductions on G = (V, E, v) then 7Z(G; Si, ... , S m ) can alternately be 
written as TZem(G; So — Si, . . . , S m _i — S m ) for V = Sq. This notation is meant to 
indicate that at the ith reduction we remove the vertices S^_i — Si from the graph 
lZem(G; So — Si, ... , S;_2 — Si_i) for each 1 < i < m where S_i — So = 0. 

With this notation in place we give a proof of theorem 13.91 
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Proof. If G = (V,E,w) G Cr and {vi, . . . ,v m } C V then, by the discussion pre- 
ceding theorem l3.8| the graph lZem(G; {vx}, . . . , {v m }) is well defined. Moreover, 
for any 1 < i < m, let G, = lZem(G; {v±}, . . . , where G\ = G. Then by 

lemma WM it follows that 1Zem[Gi; {vi}, = 1Zem(Gi; {vi + i}, fa}) which in 

turn implies that 

Hem(G; {vi}, fa}, fa +1 }, . . . {v m }) = 

Hem(G; {vi}, fa +1 }, fa}, . . . {v m }). 

By repeatedly switching the order of any two vertices as above it follows that for 
any bijection b:{vi,..., v m } ->■ {vi, . . . , v m } 

(21) Kernfa; fa},... {v m }) = Kem(G; {bfa)}, . . . {b(v m )}) . 

Suppose Si , . . . , S m induces a sequence of reductions on G. If Sq — V and each 
Sj_i — Si = {v\, . . . } vl } then theorem 13.61 implies 

Uem(G; So — Si, ... , S m -i - S m ) = 

Kem(G; fa}, {<}, . . . , {v?}, . . . , {«*}). 

If Ti , . . . , T p also induces sequence of reductions on G where Tq = V and each 
Tj_i — Ti — {v{, ...,«■} then similarly, 

TZem(G; T — Ti, . . . , T p _i - T p ) = 

^emCG;^ 1 },...,^ },...,{«?},.. .,{<}). 

If S m = T p then S m = U™ i = ULi {Ti-i-Ti). There is then a bijection 

6 : S m -> S m such that 7eem(G; {«J}, . . . , {«^}) = Kem(G; {&(«!)}, . . . , {6(«P p )}) 
implying 1Z(G; Si, ... , S m ) — TZ(G; T\, . . . , T p _i, S m ) by use of equation (12T1) . 
This completes the proof. □ 

The notation lZs m [G] = 1Z(G; Si, . . . , S m ) is then well defined for any graph 
G = (V,E,uj) in Gyr and nonempty S m C 1/. Moreover, if for any G E 0^ — 0, 
r(G) C V is nonempty then the relation G ~ i? if 7^ r (G)[G] — 7^ T (if)[-ff] is an 
equivalence relation on G w . To see this, note that as the rule t specifies a unique 
nonempty graph Tt T (G) [G] for any G G — 0, then ~ is an equivalence relation 
given that ~ is clearly reflexive, symmetric, and transitive. Hence, theorem 13.111 
holds. 

We now give a proof of theorem 14.21 

Proof. Suppose Xs{G) = (V,E,fi) is a branch expansion of the graph G = (V,E,ui). 
By assumption then each vertex of G and Xs(G) belongs to a branch of Bs{G) and 
Bs(Xs{G)) respectively. 

Suppose that e,j € E. If i = j then, as the vertex Vi belongs to some branch 
/3 G Bs{G), then u(eij) belongs to the weight sequence CIg (/?)• If i j then there 
are branches Ui, . . . , Vi, . . . , v a and u\, . . . ,v,j, . . . ,u t in Bs(G) containing Vi and 
t> j = Uj respectively. Hence, (3 = v\, . . . , v^, Vj, . . . , Ut G Bs(G). Therefore, w(ejj) 
belongs to the weight sequence Og(/3). 

Hence, each edge weight of G belongs to a weight sequence of some (3 G Bs(G). 
By similar reasoning each edge weight of X${G) belongs to a weight sequence of 
some (3 G Bs{Xs{G)). Since Bs{G) ~ Bs{Xs{G)) then uj{eij) is an edge weight of 
G if and only if uj(e{j) is an edge weight of Xg{G). Therefore, G and Xs(G) have 
the same set of edge weights. 
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Suppose p G Bs(X s (G)). As B S {G) ~ B S (X S {G)) there is a bijection 

b : B S (G) ^ B S (X S (G)) 
such that Cl Xs{G ){P) = ^gO" 1 ^))- If P(i) = v $ S define 

b(v) = (b- 1 (p))(i). 

To see that b is well defined suppose v G V — S. As u is assumed to belong to some 
branch /3 S Bs(Xs(G)) and v ^ S then u is an interior vertex of /3. Given that the 
branches of Ss(^s(G)) are pairwise independent then the branch p containing v is 
unique. Hence, b is well-defined. 

Moreover, as \P\ = \b~ 1 (P)\ then b maps interior vertices of P to interior vertices 
of b^ 1 (P). Since each vertex of V— S and V~S belong to some branch of Bs(Xs(G)) 
and Bs{G) respectively then b :V — S — >V — S \s onto. Let 

Vf = {u€V-S: b(v) =v e e V}. 

Note that \Ve\ = rit is then the number of branches in Bs(G) containing vg. More- 
over, V — S is the disjoint union 

V - s = \J v e . 

v,_£V-S 

As b : Bs{G) — > Bs(Xs(G)) preserves weight sequences then \x{e^) — uj(eu) for 
each Vj G Vi- Hence, 

(22) n ( /i (e J . i )-A)= n (n(^)- A ))= n n^)-^r- 

If S — {v\, . . . , v m } then by assumption 

~Sy(* s (G);S) for l<i,j<m. 

Equation © then implies that the edge in both lZs(G) and "^(.^(G)) have 
the same weight. Therefore, K S (G) = Us(X s (G)). As 

det (M{K S {G)) - XI) = dct (m(TZ s (X s (G))) - XI) 

then equation (l20l) implies 

det (M(G) - AJ) det (M(Af s (G)) - A/) 
det (M(G| S ) - A/) ~ det (M(X S (G)\ § ) - XI) ' 

Via equation ([7} it then follows that 

det (M(* S (G)) - A/) = det (M(G) - A/) (M^) - A) 

From equation (|22|) we then have that 



□ 



det (M(X S {G)) - XI) = det (M (G) - A/) J] (w( eii ) - A)' 

v z ev-s 



For theorem 14.51 we give the following proof. 
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Proof. Let G = (V, E, uj) and S £ sto(G). Suppose the edge weights of G are in the 
semiring U C W[A]. By theorem 14.21 both G and Xs(G) — (V,£,/i) have the same 
set of edge weights. Hence, Xs(G) has edge weights in U. 

Let ys(G) = (V, £ , v) and suppose f3 = v\, . . . , V£ £ Bs(Xs(G)), Then equation 
(jTTj) gives the weights of the edges belonging to (3 in ys{G) as 



(23) v(e 



nLi^( e fc,fc+i) for i = 1 

1 2<i<£-l 



for t ^ 1. If € = 1 then ^(en) = /i(en). Since 1 e U and U is closed under 
multiplication then the edges of f3 in ys(G) have weights in U. Given that each 
edge of ys(G) belongs to a branch f3 £ Bs(ys{G)) then this implies that ys{G) 
has edge weights in U. As Z' S {G) is a restriction of the graph 3^s(G) it follows that 
Z' S (G) also has edge weights in U. 

To construct Zg(G) from the graph Z' S (G) we add edges with weights from the 
weight set of ys(G). Hence, Zg(G) has weights in U. Suppose Zg(G) has parallel 
edges ex,...,€N of the form (a, ej) with weights wi, . . . , wn G U respectively. As 
U is closed under addition then wi + ■ ■ ■ + wn £ U. Hence, replacing ei, . . . , ejv by 
dj where eij is given the weight wi + ■ ■ ■ + wn implies that the graph Zs(G) has 
edge weights in U. 

Observe that as Bs{Xs{G)) = Bs(ys(G)) then equation (|23|) implies 

VM= II Ke k , k+ i)/\ m - 2 =V V {(3) 

k=l 

for all e Bs(X s (G)),B s (y s (G)). Hence, K S (X S (G)) = 1Z s (y s (G)). The claim 
then is that K s (y S {G)) = Ks(Z s {G)). 
To verify this let 

Bi j = {p€B lJ (y s (G);S):\(3\=e}. 

Moreover, for (3 S Bs(ys(G)) let V\,{fi) be the first edge weight of f3. Hence, if 
C = max{|/3| : (3 G B s (y S (G))} then 

(24) 2 v »w = E ( E ^w) = E (( E ^(/3))A'- 

Suppose S = {vi, ...,v m } and /3 J = vf, . . . ,vi for 1 < j < m. Following the 
^-construction we add an edge of the form (vi,Vk-\p\+2) to the graph Z' S {G) with 
weight Vl(P) for each branch in Bfj not equal to f3 J . Hence, the graph Z${G)\s 
consists of the interior vertices of /3 1 , . . . , /3 m . As these vertices have no loops and 
/3 1 , . . . , (3 m are pairwise independent then S E sta(Zs(G)). 

Suppose that Bf- — {(3i, . . . , (3n}- We then add an edge of the form (uj, Vk-e+2) 
to the graph Z' S (G) with weight Vl(f3 t ) for each branch (3 t £ Bij(yg(G)) not equal 
to ft . Hence, in Zg{G) there are N edges ei, . . . , ejv of the form (vi, vl_ e+2 ) having 
weights Vl((3i), . . . , Vl((3 m ) respectively. 

Define the branch set 

Bj j = {f3£B lJ (Zs(G);S):\f3\=l}. 
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If [} = y i) v 3 k _ e+2 , . . . , v° m then by construction Bfj — (3. Moreover, for Zs{G) — 
(V — U, £, D), the branch f3 has weight sequence 

Vz s {G)0) = E ^(/3),0,...,0,1,0,...,1 

as the weights of the edges e±, . . . , ejv are summed in the construction of Zg(G). 
Therefore, 

(25) E ^08) = E p v09). 

Note that the sum 

E t>m = E ( E w)) = E (( E ^))M*- 2 

Combining this with equations and ([25)1 implies 

E W)= E W)- 

PeBy(y S (G);S) /3GB i3 (2 s (G);S) 

Therefore, TL s iy s {G)) = K S {Z S {G)) verifying the claim. 

As 7e s (G) = K s (y s (G)) then ft s (G) = K S (Z S (G)). As S is a complete 
structural set, corollary [5] implies that G and Zs(G) have the same nonzero spec- 
trum. □ 



6. Improving Stability Estimates via Dynamical Network Expansions 

A major obstacle in understanding the dynamics of a high dimensional system is 
that the information needed to do so is spread throughout the various system com- 
ponents. By extending the notion of an isospectral branch expansion to a dynamical 
network it is possible to modify the structure of a dynamical network (J 7 , A) in a 
way that preserves the systems dynamics but concentrates this information. 

In 6.1 we introduce the concept of a dynamical network expansion (XsJ- ', Xs) of 
the network (J 7 , A). We then show that if (XsJ 7 , Xs) is globally stable then so is 
(J 7 , A) (see theorem 16. 5ft . Using this connection we show in theorem 16.71 that the 
global stability of a system (J 7 , A) is better understood by applying theorem [23] to 
(XsJ 7 , X s ) rather than (J 7 , A). 

6.1. Dynamical Network Expansions. Suppose (J 7 , X) is a dynamical network. 
Following the discussion in section 2.3, here we consider (J 7 , A) as a dynamical 
network with no local dynamics. That is, the component 

J 3 :0I^I for j el ={!,..., n}. 

ieij 

Alternatively we write 

Fj({ x i}) = Fj( x h >■•■■> x Jm)> where 1 } = {j u . . . ,j m }. 

Note that if we replace the variable Xj i of J 7 ,- by the function f(yi, . . . , yu) the result 
is the function 



3~3 ( X jl ' ■ • ■ ' X ji — 1 ' f\Vli •••)!/&)) X ji + 1 > • ■ • ? X jm ) 
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having variables Xj^, . . . , Xj i _ 1 , yi, . . . , yk, Xj i+1 , . . . , Xj m . Additionally, if the se- 
quence 7 = £1, . . . , £n let 

•^i;7 ( x ji ' • • • ' x jm ) — (xj 1 ^ , • • • , Xj m 

That is, the variables of the function are indexed by the sequences ji,£i, . ■ ■ ,£n 
for 1 < i < m. 

As defined in section 2 the graph of interactions IV = {V, E) is an unweighted 
graph. However, r> can be considered to be the weighted graph r> = (y, E, lo) G G 
by setting ui(eij) = 1 for all G E. This will allow us, in particular, to consider 
structural sets S G st(Tjr). Recall that if S — {v\, . . . , v m } then Is = {1, ... , m} 
is the index set of S. 

Definition 6.1. For S G st (Tj^) the set 

(26) A s (J r ) = {£i,...,£ N :v il ,...,v iN eB s (Tr), N>2} 

is called the set of admissible sequences of T with respect to S. 

Let (J 7 , X) be a dynamical network and suppose S G sto(Tjr). For j e 1 let 
J^i) be the function 

J-j = J-j (x,j i , . . . , Xj m ) 

in which each variable Xj e is replaced by Xj e .j if jt ^ Is- 
For i > 1 let Fqj) be the function 

•^~0>* — 1) ("^7 1 ' " ' " ' 2 '7 t ) 

in which each x 1 ^ = xe li ... ! e N is replaced by the function Tt^ if ^i ^ Is- If £\ G Is 
for each 1 < £ < t then define XsJ-j = ^(j.i-i)- 

Let ■j = £\,...,£n G As (J"). For 1 < i < |t_| = AT define the A - 2 spaces 

-X"i;7 = A^ . 

Additionally, define the functions 

XsJ~i;j 1,7) — 1,7* 

By way of notation we let 

(27) Ajv-i, 7 = A 7 , XsFn-i^ = XsJ 7 ^, and xi i7 = x^ . 

Definition 6.2. Suppose S G st (Tjr) such that each vertex of Tjr belongs to a 
branch of Bs(Tjr). Let 

jeis 7_e-4 s (.F) 
Ki<\ 2 \ 

and 

^ = (0^)®( Xi-a)- 

Ki<|7| 

The dynamical network {XsF , Xs) is called the dynamical network expansion of 
(J 7 , A) with respect to S*. 
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If 7. = £1, . . . , In and N > 1 then by construction In = j £ Is, £k ^ Is, and each 
.7-4 is a function of a;f fc _ 1 for 1 < k < N. Therefore, if £\ £ Is then the vertex 
sequence v tl ,.. . ,v iN £ B,-(I>; 5). 

If £\ Is then V£ t , . . . , is part of a branch of Bj(Tjr; S) in the following sense. 
There is a branch of the form v\ , . . . , v r -\ , , w r +i , . . . , v s £ £>j (IV; S 1 ) since is 
assumed to belong to some branch of Bj(Tjr\ S). Hence, 



which contains the vertex sequence vi 1 , . . . , vi N . 

As each branch of £>s(IV) consists of a finite sequence of vertices then for large 
enough i, l\ £ Is- As this holds for each 1 < £ < t then XsTj is well defined 
for each j £ Is- Moreover, this implies that if xg lt ,.. t g N is a variable of XsFj and 
N > 2 then v^,... ,vi N £ Bj(Tj?;S). Similarly, one can show that each XsJ-j is 
well defined for j ^ Is- 

Lemma 6.3. Suppose S £ sio(IV) and j £ Is- Then x^ t ... t i N is a variable of 
XsTj if and only if either N = 1, £% £ Is, and xi t is a variable of Tj or N > 2 
andvi 1 ,...,vg N eBj(Tjr;S). 

Proof. Suppose j £ Is and that Xi lr .. t e N is a variable of XsJ-j. A variable of XsJ-j 
is indexed by a sequence of length N = 1 if it is not replaced by a function in the 
first step of the construction of XsJ-j. Hence, Xi t is a variable of XsTj if X£ t is a 
variable of Tj and £\ £ Is- 

Assuming N = 2, or Xt lt t a is a variable of (XsT)j, then £i £ Is as cc^/a is not 
replaced by the function Tg r ^^ in the next step of the construction. However, this 
implies that x^ is a variable of J-( 2 . Hence, in constructing (XsT)j the variable 
xi t is not replaced by xt x i 2 contradicting the assumption that N = 2. Therefore, 
N 7^ 2. However, if N > 2 then as observed u^, . . ., vi N £ Bj(Tjr; S). 

Conversely, suppose j3 = v^,... ,vi N £ Bj(Tjr- S). Then xi i _ 1 is a variable of 
for 1 < i < N where £\,£n £ Is, and ^2, ■ • ■ ^ 2s- Hence, xt u ,„ t i N is a 

variable of AsJ^ x . If TV = 2 then W£ 1 , £ Bj (IV; S). As £%, 1% £ Is then xi 1 is a 
variable of Tj . □ 




if 1 Af s f 

Figure 11. The graph of interactions IV and IV s .f of (J 7 , X) and 
(AfjJ 7 , Xg) in example [T51 



To see that (XsJ-')j is well defined for each j £ Zg suppose 



X) where T is given by 
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The graph of interactions Tjr = (y,E,uj) is shown in figure [771 (left) where each 
edge is given unit weight. 

Note that S — {t"i, ^2} is a complete structural set ofTjr. Moreover, as 

Bs(I>) = {01>02> 03, &} 

where fi\ = Vi,V4,v 2 ; $2 = V\, W3, U4, «2j 03 — V\,V2] and /3 4 = 1)2, V\ then each 
vertex ofTjr belongs to a branch of Bg(Fjr). Hence, XsT is defined. 

To construct this expansion we first consider the components of J- indexed by 
Is = {1,2}- For T\ = Txix-i) note that X2 is indexed by an element of Is- Hence, 
— J~i(x2)- As each variable of Ftii) is indexed by a sequence beginning with 
an element of Is then 

X S T\ = Ti{x 2 ). 

For J-~2 = J-"2(xi,Xi) the variable x 4 is not indexed by an element of Is- Hence, 
•F(2,i) = J~2(xx,X4a)- Replacing X42 by the function J-4^2 = ^4(^142)^342) yields 
J~(2,2) = Fa{x\, J 7 4(xi42, X342)). Continuing we replace X342 by 7 r 3 ; 342 which gives 
us J~/2,3}- Since each variable of ^"(2,3) is indexed by a sequence beginning with an 
element of Is then XsT 2 = -772,3) which is given by 

XsT 2 = ^2(^1, J r 4(a;i42,-7 : 3(a;i342)))- 

Note that As (J 7 ) = {142, 1342}. For 142 £ As(J~) there is a single corresponding 
function XsJ r 2;U2 given by XsT2;\<a(x\.\/a) — xi ; i42. Similarly, the functions 

<% , S-7 r 2;1342(xi ; i342) = Xl;1342 and XsJ 7 3;l342 (X2;1342) = X2 ; 1342 

correspond to 1342 6 As (J 7 ) 

By the identification in p?7] ) the expansion (XsJ 7 , X3) is given by 



(28) X S jF{x) 



XsFi{x2) 

XsF2{Xl,X 142 ,£1342) 
Xs^U2{xi) 
^S-7 : 'l342(a;2;1342) 
^S-7 r 2;1342(xl) 



J 7 l(x2) 

?2 (Xl , Ti (X 142 , J3 (^1342 ) )) 

ari 

X2;1342 
Xl 



where Xs — \X\ ®X 2 ) © (^142 ©-^1342 ©-^2; 1342) ■ L/ie <?rap/i of interactions T x s J 7 
is shown in figure[T7\ (right). Note that the branch expansion Xs(Tjr) = Tx s f (see 
lemma\6.i 



A natural question to ask is whether the expansion (Xs T , Xs) and (J 7 , X) have 
similar dynamics. To address this let C — max{|/3| : (3 £ Bs(Tjr)}. For x £ X let 
xs £ Xs be given by 



(*s)j 



■F/(x) 

-£-i+l/-~1 



for j £ I s 



(J^- J+i (x) for j = i;h,...,£ N 
Lemma 6.4. Lei (As J 7 , Xg) be a dynamical network expansion of (J-,X). Then 

X s T*(x s )= .7f +fe (x) 
/or any x £ X, j £ X, and fe > 0. 

Proof. Let j £ Is and suppose = Tj(xj 1 , . . . ,Xj m ) is a type-1 function. Then, 
XsJ-j = J-j and jg £ Is f° r 1 < £ < m. As (xs)j e = Ff t {x) for ji £ I5 then 

(29) XsF^s) = ?i {Tf x (x), . . . , Tf m (x))) = J"/ +1 (x) 
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Suppose Tj = J 7 j{xj 1 , . . . , Xj m ) is a type-2 function. Then either ji g Is or Tup) 
is the function Tj in which the variable Xj e is replaced by the function Tj e ;j e j for 
all je Is- This implies that 

(30) (xsk=-^(*) if Jeel s . 

If on the other hand jg ^ Is then Tj e is a type-1 function. By lemma l6~3"l if 
Xi.j e j is a variable of Tj e -j t j then v^, Uj,,, i>j 6 j6j-(IV) implying i,je,j E As^). As 
per our notation, if t\, . , . ,1^ g As^) then the variable X£ lt ,„^ N — XN-i- t t lt ... t e N . 
Therefore, each variable of Tj e -j e j has the form Xij e j — ^t2;i,u,j)' Hence, 

(31) FjtUeAZs) = ^(^ £ - 2+1 (x)) = j£(x) for u i Is- 

As the variable xg . is replaced by the function Tj e [ j e j for all ji ^ Is then (|30j) 
and (|3"Tj) together imply that 

AsT) (x s ) = 7} (j* (J) , . . . , J-/ m (x)) = (x) 

as in (|29|) . Continuing in this manner it follows that 



(32) XsFj(is) = X s Tj{St) = Tf +1 (x) for j g Z s . 

Hence, lemma [6~4l holds for fc = 1. 
Note that each 

XsTi^is) = X s T i .,i 1 ,..., iN (Tt ii - 1)+1 (*)) = rfr +2 W- 
This together with (13"2l implies that 



^ +i (x) forigls 



^" l+2 (x) for j = i;£ 1> ...,£ 



N 



As XsT(x.s) has the same form as xs then replacing xs by XsTfes) and ap- 
plying the same argument yields 



^(x s ) 



J-/ +1 (x) for 3 els 
j£- l+2 (x) for J=i;*i,...,4v 



Continuing in this manner it follows that X s T^(5ts) = -7\f +fc (x) for all j g X 5 , and 
k > 0. □ 

Theorem 6.5. Suppose {XsT , Xs) is a dynamical network expansion 0/ (J 7 , X). 
If (XsT, Xs) has a globally attracting fixed point then (T, X) also has a globally 
attracting fixed point. 

Proof. Suppose y g Xs is a globally attracting fixed point of the dynamical network 
(X S T,X S ). For x g X and j g I s lemma implies that X s Tj(x) = jf +fe (x) 
for all fc > 0. Hence, for j g I5 

(33) d(Tf(x),yj) ->-0 as fc -> 00. 

Proceeding as in the proof of lemma 15741 suppose j ^ Is and that the function Tj 
is a type-1 function. Hence, Tj = J-j (x^ , . . . , Xj m ) where ji g Is for all 1 < I < m. 
Let e > 0. As the function J 7 : X — >• X is continuous and 

Jj £ (x) = 7- i (j*- 1 (x),...,7^- 1 (x)) 
then (j3"3")l implies, for large enough k, that d(Tj(5c),Tj(y)) < e. Therefore, 

(34) d(T^(iL),Tj(y)) -)■ as fc -4 00 
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if J-j is a type-1 function. 

Suppose j £ Ig and that the function J-j is a type-2 function. If variable 
of J-j then either jj> G Is or J 7 ^ is a type-1 function. If jg G Is then (|33f> implies 
that J 7 * (x) — > y Jf as — > 00. If je. ^ Is then J 7 ^ is a type-1 function and 
j£(x) -> J^(y) by {35). Hence, J^x) -> JJ°(x) elasfc^oo where the 
component J"?°(x) depends only on y. 

As each J 7 , for j G I is a type-i function for some i < 00 then continuing in this 
manner implies that (J 7 , X) has a globally attracting fixed point. □ 

The following is an immediate corollary to theorem 16.51 

Corollary 6. Let (XgJ-, Xg) be a dynamical network expansion of (J 7 , X) . Suppose 
the constants Ajj satisfy (0) for As J 7 . // /f(A) < 1 then (J 7 , X) has a globally 
attracting fixed point. 

The question then is whether there is any advantage in considering a dynamical 
network expansion over the original (unexpanded) network. As it turns out, dy- 
namical network expansions always allow for better estimates (or at least no worse) 
of a dynamical network's global stability. To demonstrate this we first require the 
following lemma. 

Lemma 6.6. Suppose that (As J 7 , As) is a dynamical network expansion o/(J-", X). 
Then X s (Tj?) = T XsT . 

Proof. Let Tx s F — (V, By construction 

V = S U K 2 : 7 6 As(I), K i < 111} 

where the vertex Vj G S corresponds to the component XgTj and Vi- n G V — S the 
component XgTi a . 

Suppose 7 = ti, . . . , £ N G As (J"). Then v tl ,...,vt N G S,s(I» where A > 2. 
Hence, lemma 15731 implies that xt lt ... t i N — £at_i i7 is a variable of the function 
XgTg N . Moreover, as XgTi- ri = XgJ 7 i n (xi-i n ) for 1 < i < 7 where xi n = Xi t then 
the vertices vg 1 , V2--y, • ■ • , Vjv-1,7, form a path in TxsF such that , vi N G S". 

As any two admissible sequences of Ag(F) are distinct then Tx s r\g consists 
of |As(.F)| paths that share no vertices of the form «2 l7 , . . . , vn-i,j- Since the 
vertex v in G V has no loop in ^XsMg then Tx s f\§ contains no cycles. Hence, 
S G sto(Tx s r) and for each 7 = t\, . . . ,£n, 

(35) Bi 1: £ N {T Xs jr; S) = {v ei ,v 2 , 2 , ■ ■ .,v N -i a ,V£ N : 7 = ii,. ..,l N G Ag{T)} 
if A > 2. 

For /? = Vt 1 ,...,vt N G B.e lr e w (I>; S) and j = h,---jN let 
'/3 if A = 2 



otherwise 



(36) b(J3) 
The claim is that 

b : Be 1 ,e N {Tj?; S) -> S^^„(r^ g jr; 5) 
and that this map is a bijection. 

To see this suppose A = 2. Then, (3 = V£ 1 , G Bi lt £ N (T^; S) implying xi 1 is a 
variable of Tg 2 . As i\ G Is then by construction x^j is a variable of XgTg 2 . Since 
u^D^a G S then /3 = v tl ,ve 2 G £>^,£ 2 (r^ s jF; 5). Conversely, suppose v il ,ve 2 G 
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Be 1 ,e 2 {^'x s J r ] S) i.e. xi x is a variable of the function XsJ-"e 2 . Then X£ t is a variable 
of Te 2 by lemma HT3l implying 8 G Bi 1 ^ 2 (Tjr; S). Hence, b restricted to branches of 
length 2 is a bijection. 

For N > 2, the sequence 7 = l\, . . . G -As^-F). Hence, 

6(/3) = V tl , U2;7; • ■ • , WJV-1;7> V£ N 

which is a branch in Bi lt i N (r^ s jr; 5) by (|55j) . Conversely, if 

%j "2;7, • ■ • , vn-i-j, V£ n G ^^^(r^jr; S) 

then 7 £ As^) implying vi 1 ,...,vt N G B£ 1 .£ N (Tjr; S). Therefore, the function 
b : Bi l .i N (Yjr;S) -)■ Bt u t N (r^ s ^; S) is a bijection. 

As £i,£n G Is were arbitrary and each edge of IV an d r^ s jr have unit weight 
then 

B S (I» * B s (T Xs r). 
Since the branches of r ' x s T are independent it follows that Xs(Tjr) = Y x s f- □ 

Theorem 6.7. (Improved Stability Estimates for Dynamical Networks) 

Let {XgJ-, Xg) be a dynamical network expansion of (J-,X). Suppose there exist 
constants Ay satisfying 0] for J 7 . Then there are constants Ay satisfying (0) for 
X S T such that p(A) < p(A). 

Before proving theorem 16.71 suppose the constants Ay satisfy @ for the map 
T. Define Tjr(A) = (V,E,ui) to be the graph with adjacency matrix A where Ay- 
satisfy © for J 7 . Note that as Ay ^ if and only if Fj is a function of Xi then 
Tjr and Tjr(A) are identical as unweighted graphs. Therefore, if S G sio(IV) then 
S G ato(I>(A)) and B s (l>(A)) = B S (I». 

Moreover, the branch 8 = Wy, . . . ,Vi m G Ss(IV(A)) has weight sequence 

^rv(A)C#) = Ayy, . . . , 0, A i;j ^ +1) 0, . . . , A im _ 1 , im 
implying that the branch product 

(37) P u (fi)= [[ %^T +1 - 

We now give a proof of theorem 16.71 

Proof. Let xjel s . For 7 G As {J 7 ) and 1 <i <\j\, 

d(X s J'i a (^), Xs^jiy)) = d(xi-i- :L ,yi-i- 1 ). 
Therefore, for any 7 G As (J 7 ) and 1 < i < I7I the constants 



(38) A 



if j = 2; 7 
otherwise 



satisfy equation ([2]) for Xs J-j . 

For j £ I let if = Ij U I5 and IJ = Ij — Is- For ii G Xg 

(39) d^sJ^xJ.^fy))^ 

^ ^ -^-22*1 d{xi 2 , y 22 ) -H ^ ^ K.j i ^j iX d{j~i 2 - i) j iX (x), ^*i 2 ;ii (y)) * 
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Similarly, the sum 

(40) A l2ll 4j- 42;ll (x), j; 2;n (y)) < 

^ ' Aj 2 ii f ^ ] A-t3»2^('^*3,»2,in 2/«3,»2,ii ) + ^ ] Aj 3 j 2 ^(•^i'sife ,»i ( x ) ) "^"isjta ,» i (y)) 



Note that if i<x G 1^ and i% G 2j then ii ^ Is and 13 G Is implying the sequence 
v ia ,v i2 ,v ll G Bs(l>(A)). In particular, ^(v^ , w l2 , t^JA = A l3 j 2 A l2 , i:1 by 

Let ={/3e B^ (I>(A); 5) : \/3\ = £}. Moreover, if the branch /3 = v 1 ...,v m 
let xp = x\ m and 7(/3) = l,...,m. If Tjr(A) = (V,E, uS) then combining in- 
equalities in (|39| and (|40|) yields 

d(-Y s J' il (x),A'sJ- il (y)) < 
' Aj 2 jj d{xi 2 , j/i 2 ) + S ' S ' Ai 3 j 2 A{ 2 i 1 d(d(xi 3i j 2! j 1 , j/j 3) i 2) i x ))+ 

^ ' ^ , Ai 3 i 2 Ai 2 i! d(j r i 3; i 2j i 1 (x), ^ r i 3 ;i 2 ,ii (y)) = 

)+ E (^(/3)A)d(^, y/3 )+ 

^ ] ^ ] Ai 3 i 2 Ai 2 ij d(^J-i a] i 2t i 1 (x), J~i 3 -i 2 ,i 1 (y)) • 



Continuing in this manner it follows that 

c 

]T ^ 2 * 1 d{x l2 ,y %2 ) + Y J ( (W)^ l ~ 2 )d(x ,yp)) = 
^ ' Ai 2 i! d(xi 2 , j/i 2 

)+ E (^(/3)A I/31 2 )d{x ,y ) 
i 2 ex+ peB^iXA)^) 

where £ = max{|/3| : /3 G B<?(IV(A))}. 

By lemma 16.31 the variables of X$ are indexed by either j G Is or by the 
sequence j(/3) for some (3 G Bs(IV). Hence, 

{Aij if i, j G l s 

V U (P)X^- 2 if j G Is and i = 7 (/J) 
otherwise 

satisfy condition ([2]) for XsTj. The constants in (|38|) and (|41 [) therefore give a 
complete set of the constants A^- satisfying ^ for As-F- 

Suppose S — {v\, . . . ,v m }. Note that the bijection defined in ([31)1 also gives a 
bijection 

(42) 6 : Bij(Tj^(A); S) Bij(T Xs ^(A); S) for l<i,j<m. 
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In particular, if (3 € i3s(IY(A)) then equations (|38| and (|4Tj) imply that 

°r^(A) W)) = 1, 0, 1, . . . , 0, V^P)\^- 2 . 
Let r^ s ^(A) = (V,£,(J,). As \/3\ = \b(/3)\ then 

Given the one-to-one correspondence noted in equation (|4"2")l it follows that 

As 5 <E st (I>(A)) and 5 G si (rA , s jf(A)) then corollary [3] implies that I>(A) and 
r^ s ^-(A) have the same nonzero spectrum. Hence, p(A) = p(A) for this choice of 
constants completing the proof. □ 



Example 14. Let C{x) = Ax(l— x) be the standard logistic map and let the junction 
Q(x) = l~x 2 
systems 



where both C : [0, 1] — > [0, 1] and Q : [0, 1] — > [0, 1]. Suppose the local 
Ti{xi) 



C{xi) for i = 2,4 
Q(xi) for i = 1,3 



Let F : [0, l] 4 



[0, l] 4 be the interaction given by 
l , 



F(x) 



1^4 



^2 

4^1 
4^1 + 



Consider the dynamical network (J 7 , X) where T = F oT and X = [0, 1] . Note 
that the constants Li = 4 for i = 2,4 and Li = 2 for i = 1,3 satisfy equation (QJ) 
/or t/ie ZocaZ systems (T,X). For Ay = max xe x | (DF)ji (x) | £/ie matrix 



i i 

4 4 





4 



_L 

4 

1 




1 

4 



diag[2,4,2 7 4] 



1 



l 

2 





1 
10 



As the spectral radius p(Mjr) = 1.08 > 1 then theorem \2. 51 does not directly apply 
to the dynamical network T . 

Moreover, suppose J- = (F o T) o id generates the dynamical network (J 7 , X) 
having no local dynamics. For A^ ~ max xg x \ {DF)ji(x)\ one can compute that 
A = Mjr. As A > 1 then considering (J-,X) as a network without local dynamics 
does not improve our estimate of the network 's stability as it did in example [U 

However, note that 



F(x) 



Fi(x 2 ) 




T\ {x\,x£) 




■Fi(xi) 




_ T\{x\,xz) 





|G(*3) 
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As T has the form of the map given in example UM then following this example 
one can compute 



XsJ-i 





















\C(x 2 ) 



\Q{xi) + \C{Q(x li2 ) + \Q{\Q(x 13i2 ))) 



4 
Xx 

£2;1342 
Xi 



for S = {vi_, v 2 }. Here, X 142 , X 2 . 1342 , X 1342 = [0, 1] implying X s 
Letting Kij = max l6 x s | {DXs^ji (x) | it follows that 



[0,1] E 



A 



1/2 110 

1 
0.265 
1 
0.012 



As p(A) — 0.90 < 1 corollary [5| implies that the dynamical network (J-,X) does in 
fact have a globally attracting fixed point. 

Importantly, the method of using network expansions generalizes the method 
used in [1] [2] for determining whether a network has a unique global attractor. As 
a demonstration we give the following examples. 



Example 15. For the dynamical network (J 7 , X) let 



L = max Li and A r , 

Ki<n 



max 



A, 



Recall the spectral radius of the graph Tjr (considered as an unweighted graph) 
is p(Tjr}. The condition in pQ guaranteeing the global stability of the dynamical 
network {!F,X) is the following: Lf the product 



(43) 



LA 

ina:r P(IV) < l 



then (J-, X) has a globally attracting fixed point. 

Consider the dynamical network A") given in example® In this case L = 4, 
A max = 2/9 7 and p(Tjr) = 2. Hence, 

LA maxP (T T ) = 16/9 > 1. 



Therefore, (J-,X) does not satisfy condition C^l ). However, as the spectral radius 
p(A4jr) < 1 then the dynamical network (J 7 , X) has a global attractor by theorem 
731 (see example 2). 



Example 16. Suppose (J 7 , A") is the dynamical network given in example \H\ Using 
the mathematical framework established in [5] one can calculate that the topological 
pressure p( fe ) (ip) of the system is given by 

P^(ip) = ln(1.08). 

The condition given in 2 that guarantees the global stability of the dynamical net- 
work (J-,X) is if the systems topological pressure 

p [k \v) < o 
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(see theorem 2, [2],). As ln(1.08) > then the method introduced in [2] does not 
imply the global stability of (J-,X). However, by expanding the network it is pos- 
sible to demonstrate that the system does have a globally attracting fixed point (see 
example \lJ$ . 

As a final observation in this section, we note that it is possible to sequentially 
expand a dynamical network and thereby improve ones estimate of whether the 
original unexpanded system has a globally attracting fixed point. 

7. Concluding Remarks 

The goal of this paper is twofold; introduce isospectral graph transformations 
and use these transformations to investigate the dynamic properties of dynami- 
cal networks with arbitrary topologies (graph of interactions). For the first, we 
described the general process of isospectral graph reductions in which a graph is 
collapsed around a specific set of vertices. We then considered sequences of isospec- 
tral reductions and showed that under mild assumptions a typical graph could be 
uniquely reduced to a graph on any subset of its vertex set. 

In [10] , such reductions and sequences of reductions are used to improve the clas- 
sical eigenvalue estimates found in [TTJ HH [27] . These estimates of Gershgorin 
et al. can be traced back to work done by Levy, Desplanques, Minkowski, and 
Hadamard [22j [T3l E31 [17] on diagonally dominant matrices. The main result of 
[TO] is that such eigenvalue estimates improve as the graph associated to a matrix 
is reduced via the method of isospectral reduction given in this paper. 

Isospectral reductions are also relevant to the theory of networks as they provide 
a flexible means of coarse graining networks (i.e. viewing the network at some ap- 
propriate scale) while simultaneously introducing new equivalence relations on the 
space of all networks. From the point of view of applications, isospectral graph re- 
ductions are a tool which allow experimentalists to modify the networks with which 
they are working. However, the particular transformation requires the expertise of 
the experimentalist (biologist, ecologist, etc.) to determine an appropriate set of 
network elements (vertices) over which to reduce the network. 

In this paper we additionally considered isospectral graph transformations over 
fixed weight sets. Such transformations allow one to transform an arbitrary graph 
while preserving the set of edge weights of the graph along with its spectrum. Mo- 
tivated by these procedures we also introduced the notion of a dynamical network 
expansion. 

Dynamical network expansions modify a dynamical network (i.e. a dynamical 
system with a graph structure) in a way that preserves the system's dynamics 
but alters its associated graph structure (graph of interactions). In this paper 
we demonstrated that this procedure allows one to establish global stability of a 
more general class of dynamical networks than considered in [TJ [2]. That is, such 
expansions provide more information and can be used to obtain stronger results 
regarding network dynamics than the direct analysis of the network. It is worth 
mentioning that, as any finite dimension dynamical system with discrete time can 
be considered a dynamical network, this technique can be applied to a very large 
class of systems. 

Lastly, the results of the present paper introduce various approaches to sim- 
plifying a graph's structure while maintaining the graph's (network's) spectrum. 
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Therefore, these techniques can be used for optimal design, in the sense of struc- 
ture simplicity of dynamical networks with prescribed dynamical properties ranging 
from synchronizability to chaoticity [TJ 0] ■ 
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